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I.  INTRODUCTION 


This  chapter  addresses  the  followirtg  issues; 

1 .  The  motivation  behind  the  research. 

2.  The  objectives  of  the  research. 

3.  A  preview  of  subsequent  chapters. 

A.  MOTIVATION 

Since  1982  the  Naval  Postgraduate  School  (N'PS)  has  been  examining  in  depth  the 
inventorv^  control  models  in  use  at  the  Navw's  Inventory  Control  Points  (ICP).  The  two 
main  ICP's  are  the  Ships  Parts  Control  Center  (SPCC)  and  the  Aviation  Supply  Ofiice 
(ASO). 

This  work  was  undertaken  at  the  request  of  the  Naval  Supply  Systems  Command 
(NAVSL'P)  in  conjuction  with  NAVSL'P's  Resolicitation/ Resystemization  Project.  The 
primary  objective  of  the  resolicitation  project  was  the  acquisition  of  new  computer 
hardware  for  the  ICP's,  but  the  effort  also  provided  the  opportunity  to  make  software 
improvements  that  would  make  better  use  of  the  new  hardware. 

Initially,  NAVSUP  asked  NTS  to  "develop  improvements  to  the  existing  peactime 
wholesale  provisioning  models  for  secondarv^  items  used  by  ...  SPCC  and  ...  ASO."  [Ref. 
1:  pp.  1]  Provisioning  models  are  those  that  determine  the  range  and  depth  of  repairable 
and  consumable  items  that  should  be  initially  procured  in  support  of  a  new  weapon 
system.  A  provisioning  model's  objective  is  to  procure  enough  spares  to  provide  ade¬ 
quate  support  for  the  weapon  system  until  the  first  w'holesale  replenishment  order  ar¬ 
rives. 

In  response  to  NAVSUP's  initial  request,  NTS  proposed  several  provisioning  models 
(Ref.  1]  and  conducted  detailed  comparisons  between  these  models  and  those  being  used 
by  the  ICPs  at  that  time.  (Ref  2]  As  a  result  of  these  analyses  a  new  provisioning  model 
was  adopted  at  the  ICP's  in  1984.  This  model  is  called  the  Mean  Supply  Response  Time 
(.MSRT)  model  because  its  objective  is  to  minimize  the  MSRT  for  a  set  of  secondary 
items  belonging  to  a  new  weapon  system.  Mean  supply  response  time  is  defined  as  the 
mean  time  it  takes  the  supply  system  to  satisfy  the  demand  for  an  item. 

The  logical  next  step  was  to  develop  a  replenishment  model  having  the  same  objec¬ 
tive  function.  This  model  would  be  used  to  make  all  subsequent  buys  of  inventor)’  after 
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the  initial  provisioning  buy.  The  initial  effort  in  this  direction  concentrated  on  a  re¬ 
pairable  replenishment  model. 

The  initial  attempt  to  extend  the  MSRT  concept  to  repairable  item  management 
was  described  by  Gormly  in  his  master's  thesis  at  NFS.  [Ref.  3  )  His  thesis  investigates 
an  aggregate  demand  inventory  model  for  repairable  items  that  assumes; 

the  probability  distribution  for  the  inventorv’  position  (defined  as  on-hand  plus  on- 
order  plus  in-repair  minus  backorders)  was  ...  uniform  with  its  equally  likely  states 
being  a  function  of  a  weighted  sum  of  the  procurement  quantity  Qp  and  the  repair 
batch  induction  quantity  Qi.  [Ref  4;  pp.  3j 

This  average  was  then  incorporated  into  the  continuous  review  model  of  Hadley  and 
Whitin,  [Ref  4;  chap.  4).  A  different  approach  was  taken  by  Apple,  also  for  his  master's 
thesis  at  NFS.  [Ref  5]  The  model  presented  by  Apple  took  a  queuing  model  approach 
towards  evaluating  the  delay  created  by  waiting  for  carcasses  to  accumulate  for  repair, 
and  attritions  (not  repairable)  to  "accumulate"  towards  a  procurement  quantity.  The 
model  was  actually  a  provisioning  rather  than  a  replenishment  model. 

In  the  spring  of  1987,  McMasters  attempted  to  prove  the  inventorv-  position  prob¬ 
ability  distribution  assumed  in  Gormly's  thesis.  He  discovered  that  the  probability  dis¬ 
tribution  assumption  was  incorrect.  He  also  found  that  the  P„,  repairable  model 
developed  by  Fleet  Material  Support  Office  (FMSO)  had  been  based  in  the  same  as¬ 
sumption  and  hence  was  also  incorrect.  McMasters  was  able  to  develop  a  correct  ver¬ 
sion  of  the  probability  distribution  for  the  case  of  a  repair  survival  rate  of  100‘*'o  (all 
carcasses  inducted  for  repair  will  be  successfully  repaired).  Research  continues  on  the 
case  of  a  repair  survival  rate  less  than  100%. 

Fortunately,  the  problems  with  the  probability  distribution  of  the  inventoiy  position 
do  not  exist  if  only  consumable  items  are  considered.  The  inventory  position  distrib¬ 
ution  is  merely  that  of  the  exact  formulation  from  Hadley  and  Whitin.  [Ref  4:  chap  4] 
The  consumable  version  of  the  .MSRT  model  is  also  a  limiting  case  of  the  repairable 
model.  As  a  consequence,  it  is  convenient  to  conduct  a  comparison  between  the  con¬ 
sumable  version  of  the  MSRT  model  and  the  current  L'lCF  inventoiy  model  rather  than 
waiting  until  the  MSRT  repairable  model  is  completely  developed.  Such  a  comparison 
will  provide  insights  into  the  benefits  of  the  MSRT  model  and  possibly  result  in  its  im¬ 
plementation  for  consumables  rather  than  waiting  for  the  repairable  version. 
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B.  THESIS  OBJECTIVES 

1.  Develop  a  computer  program  for  the  current  ICP  model  for  the  inventor>’  man¬ 
agement  of  consumable  items. 

2.  Derive  an  MSRT  model  for  consumable  items  and  develop  a  computer  program  for 
it. 

3.  Using  ICP  provided  inventorv'  management  data,  make  a  comparative  evaluation 
of  the  performances  of  the  proposed  MSRT  model  and  the  current  ICP  model. 

C.  PREVIEW  OF  SUBSEQUENT  CHAPTERS 

In  chapter  II  a  discussion  of  the  Nav\-  Supply  System  is  provided  to  give  the  reader 
an  understanding  of  the  environment  in  which  the  models  are  designed  to  operate.  This 
is  followed  by  an  in-depth  look  at  the  current  ICP  model,  including  the  mathematical 
formulas  and  the  constraints  that  are  applied.  In  chapter  III  the  MSRT  model  is  pre¬ 
sented  and  the  algorithim  used  for  optimization  is  described.  Chapter  IV  discusses  the 
data,  hardware  and  software  used  in  the  research.  Chapter  V  presents  the  measures  of 
effectiveness  used  to  assess  model  performance,  the  results  of  the  comparative  analysis 
and  an  evaluation  of  these  results.  Chapter  VI  provides  a  brief  summar\’  and  recomm¬ 
endations. 
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II.  THE  NAVY  SUPPLY  SYSTEM  AND  THE  UICP  CONSUMABLE 

INVENTORY  MODEL 

The  intent  of  this  chapter  is  to  introduce  the  reader  to  the  supply  system  in  which 
the  two  inventory  models  that  form  the  crux  of  this  research  are  designed  to  operate. 
The  chapter  also  explains  how  the  LTCP  model  functions  within  that  system.  This  will 
be  accomplished  by: 

1.  Giving  an  overview  of  the  Na\y  Supply  System. 

2.  Describing  the  wholesale  consumable  inventory  system. 

3.  Explaining  basic  inventory  theory. 

4.  Deriving  the  equations  used  in  the  LTCP  model. 

5.  Discussing  the  constraints  placed  on  the  LTCP  model. 

A.  NAVY  SUPPLY  SYSTEM 

The  purpose  of  this  section  is  to  identify  the  structure  to  which  the  Navy  Supply 
System  belongs  and  the  components  which  make  up  the  Navy  Supply  System. 

The  Naval  Supply  Systems  Command  (NAVSLP)  is  responsible  for  the  manage¬ 
ment  of  the  Navy  Supply  System.  The  Commander  of  NAVSLP  gives  policy  direction 
and  guidance  on  inventory  management  of  all  material  within  the  system.  The  actual 
wholesale  supply  system  management  functions  are  carried  out  by  inventory  control 
points  (ICPs).  The  major  Navy  ICPs  are  the  Aviation  Supply  Office  (ASO)  for  aviation 
material  and  the  Ships  Parts  Control  Center  (SPCC)  for  surface  ship  and  submarine 
items. 

The  ICPs  provide  two  broad  categories  of  support:  program  support  and  supply 
support.  While  both  categories  are  equally  vital  to  the  success  of  the  Navy  Supply 
System,  the  supply  support  category  is  of  most  importance  to  this  thesis.  Supply  sup¬ 
port  functions  include: 

1.  Budget  development  for  parts  support. 

2.  Requirements  determination. 

3.  Material  procurement. 

4.  Integrated  inventory  management. 

5.  Material  distribution  and  issue. 

6.  Repairables  management. 
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This  thesis  will  concentrate  on  the  functions  of  requirements  determination  and  budget 
development  for  parts  support. 

The  Nav7  stock  points  are  the  holders  of  the  inventories  in  the  Navy  Supply  System. 
"The  main  mission  of  a  stock  point  is  the  physical  distribution  of  material."  [Ref.  6:  pp. 
26]  Their  functions  include: 

1.  Receiving  and  stowing  material. 

2.  Issuing  and  shipping  material. 

3.  Reporting  receipts  and  issues  to  each  item's  ICP. 

Stock  points  are  located  throughout  the  United  States  and  overseas,  near  major  con¬ 
centrations  of  operational  and  industrial  customers. 

The  two  types  of  material  managed  by  the  Navy  Supply  System  are  repairables  and 
consumables.  "A  repairable  is  an  item  of  supply  that  can  be  made  to  function  by  a  re¬ 
pair  process  after  it  breaks."  (Ref.  6:  pp.  4]  Examples  of  repairables  include  gear  boxes 
and  circuit  boards.  A  consumable  is  an  item  which,  when  broken,  during  use,  is  not  re¬ 
paired.  Examples  of  consumables  are  valves,  steam  pipe  and  flanges.  The  inventory 
models  involved  in  this  thesis  research  deal  only  with  consumables. 

The  Na\y  groups  similar  items  for  administrative  and  management  purposes  into 
Cognizance  Groups  (COG),  which  are  designated  by  two  alphanumeric  symbols  or 
"digits".  The  COG  is  used  to  identify  the  ICP  or  agency  responsible  for  managing  the 
item.  It  also  indicates  whether  the  item  is  a  repairable  or  consumable.  The  consumable 
items  used  in  this  research  have  a  "COG"  of  IH,  indicating  that  the  items  are  consum¬ 
ables  managed  by  SPCC. 

To  foster  inventory  management  and  flexibility,  COGS  are  further  divided  by  adding 
two  more  digits  to  them.  These  last  two  digits  segregate  items  by  essentiality  or  weapon 
system  and  by  requisition  frequency.  For  example,  the  4-digit  cog,  1H4A  represents  a 
IH  cog  item  having  an  essentiality  of  4  (the  highest)  and  a  requisition  frequency  of  3  or 
more  per  quarter  (denoted  by  the  A). 

B.  WHOLESALE  CONSUMABLE  INVENTORY  SYSTEM 

Within  the  Na\y  Supply  System  there  are  three  levels  of  inventorj’  management  of 
items: 

Wholesale  -  The  inventory  is  managed  by  ICP  inventory  managers.  The  inventory' 
managers  have  worldwide  visibility  and  control  of  this  inventory. 
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Retail  Intermediate  -  This  level  of  inventory  is  designated  to  support  a  specific  ge¬ 
ographic  area  and  is  normally  managed  by  the  Nav>’  stock  point  in  that 
geographic  area. 

Retail  Consumer  -  This  is  a  level  of  inventory  held  for  an  activity's  own  use. 

The  inventory  models  in  this  research  deal  only  with  the  wholesale  inventory  level. 

The  general  characteristics  of  wholesale  level  inventories  are: 

1.  Inventory  levels  are  computed  based  on  system- wide  demand  data. 

2.  The  ICP  inventory  manager  knows  where  the  material  is  located  and  has  unre¬ 
stricted  access  to  it. 

3.  The  ICP  must  approve  all  requisitions  for  items  from  these  inventories. 

4.  Material  is  assigned  to  one  or  more  retail  stock  points  for  storage  and  issue  by  the 
ICP  inventoiy’  manager. 

As  stated  in  item  1  above,  wholesale  inventory  levels  are  computed  based  on 
system-wide  demand  data.  How'  does  this  demand  data  reach  the  ICP  inventory  man¬ 
ager?  It  is  the  responsibility  of  the  stockpoints  to  submit  the  demand  data  to  the  ICPs. 
When  a  stock  point  issues  a  consumable  to  an  end-use  customer,  the  stock  point  reports 
the  issue  to  the  cognizant  ICP  via  a  transaction  item  report  (TIR).  Stock  points  also 
use  TIRs  to  report  receipt  and  redistribution  of  material.  The  information  from  the 
TIRs  is  immediately  processed  into  the  ICP  data  files.  This  and  other  information  in  the 
files  aids  the  ICP  inventory'  manager  in  making  decisions  involving  requirements  deter¬ 
mination,  material  distribution  and  procurement  of  replenishment  stock.  Data  from 
these  files  was  used  in  the  comparative  evaluation  of  the  current  UICP  and  proposed 
MSRT  inventory-  models. 

C.  THE  BASIS  FOR  THE  ICP  INVENTORY  MODEL 

There  are  two  major  questions  an  ICP  inventory  manager  must  answer  with  regard 
to  consumable  items; 

1.  How’  much  to  order;  i.e.,  the  reorder  quantity  (Q)? 

2.  When  to  order;  i.e.,  the  reorder  point  (R)? 

In  an  unconstrained  world  the  ICP  inventory  manager  could  order  as  much  as  the  cus¬ 
tomer  needed  and  order  as  frequently  as  demands  occurred.  Unfortunately,  the  real 
w'orld  is  full  of  constraints  with  w'hich  the  inventory  manager  must  deal.  These  include 
inventory  costs. 

The  variable  inventory  costs  which  are  dependent  on  Q  and  R  are  generally  divided 
into  three  categories: 
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Order  Costs  -  costs  accrued  when  determining  and  processing  buys,  and  the 
stockpoints'  receipt  and  stowage  costs. 

Holding  Costs  -  includes  costs  of  investment,  storage,  obsolescence  and  pilferage. 

Shortage  Costs  -  costs  of  not  having  items  when  needed. 

These  three  costs  plus  the  cost  of  buying  the  items  (which  is  not  dependent  on  Q  and 
R)  make  up  the  total  inventory  costs. 

Two  basic  types  of  inventory  models  exist  which  consider  these  inventory  costs  when 
computing  reorder  quantity,  Q  and  reorder  point,  R.  They  are  the  continuous  review 
and  periodic  review  models.  The  consumable  model  used  by  the  Navy  ICPs  for  whole¬ 
sale  inventory  is  a  type  of  continuous  review  model  called  the  min-max  model.  The 
min-max  model  seeks  the  values  of  Q  and  R  which  minimize  the  average  annual  variable 
inventory  costs. 

In  the  min-max  model,  a  replenishment  order  occurs  when  an  item's  inventory  po¬ 
sition  (on- hand  plus  on-order  minus  backorders)  hits  or  falls  below  R.  At  the  time  of 
ordering,  the  amount  ordered  is  Q  plus  the  difTerence  between  the  reorder  point  value 
and  the  inventory  position  value  at  that  time.  In  the  simplest  form  of  the  min-max 
model,  Q  is  usually  the  amount  which  minimizes  the  average  annual  ordering  and  hold¬ 
ing  costs.  For  consumables,  this  is  the  Wilson-Harris  economic  order  quantity  [Ref  7]: 


!%ad 
V  IC  ’ 


where; 

A  =  administrative  order  cost; 

D  =  forecasted  average  quarterly  demand; 

I  =  holding  cost  rate;  and 
C  =  cost  of  one  unit  of  the  item. 

The  reorder  point,  R,  for  the  simplest  form  of  the  min-max  model  is  determined  from 
minimizing  the  costs  of  carrying  safety  stock  and  incurring  backorders  (shortages).  It 
is  a  function  of  leadtime  demand  and  the  variability  of  demand.  The  formula  for  R  is: 


R^{DxL)  +  SL\ 
where: 

L  =  the  forecasted  mean  procurement  lead  time  in  quarters;  and 
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SL  =  safety  level,  a  function  of  demand  and  lead  time  variability. 

The  reorder  point  is  computed  by  first  computing  the  optimal  risk  of  a  stockout 
according  to  the  following  formula  [Ref  7;  pp.  145] ; 

QIC 

e-c+-4;.o  ^ 

w’here: 

=  shortage  cost  of  one  demand  backordered  for  one  year. 

The  optimal  risk  is  used  with  the  probability  distribution  for  demand  during  procure¬ 
ment  leadtime  to  derive  the  actual  reorder  point  value.  This  reorder  point  value  is  the 
mean  leadtime  demand  (Dx  L  which  is  designated  by  ii),  plus  the  safety  level.  For  the 
normal  distribution,  the  safety  level  is 

SL  =  zx  a\ 

where: 

z  =  the  standard  normal  deviate  for  the  optimal  risk;  and 
a  -  the  standard  deviation  of  leadtime  demand. 

How  else  is  the  performance  of  an  inventory  model  judged  besides  by  the  minimi¬ 
zation  of  average  annual  variable  costs?  Typically,  several  other  Measures  of  Effective¬ 
ness  (MOEs)  are  used.  A  MOE  is  a  function  of  the  decision  variables  (in  this  case  Q 
and  R)  and  inventory  system  parameters.  The  tw'o  major  MOEs  the  Navy  Supply  Sys¬ 
tem  uses  are: 

1.  How  well  the  available  Naw  Stock  Fund  (NSF)  budget  is  allocated  (another  cost 
type  MOE). 

2.  Supply  material  availability  (SMAXcustomer  service  type  of  MOE). 

In  this  era  of  diminishing  resource  dollars,  the  Navy  Supply  System  must  allocate 
its  NSF  inventory  investment  dollars  to  maximize  customer  support.  The  UICP  model 
does  this  by  allocating  these  investment  dollars  to  the  optimal  order  quantities  and  re¬ 
order  point  which  minimize  total  average  annual  inventory  costs  while  attempting  to 
achieve  a  minimum  goal  of  85%  SMA.  The  dollar  investment  needed  to  support  the 
inventor}'  depth  (order  quantity  +  reorder  point)  is  used  by  NAVSL’P  in  formulating 
NSF  budget  requirements. 
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D.  THE  UICP  WHOLESALE  CONSUMABLE  INVENTORY  MODEL 

In  this  section,  the  inventory  model  used  by  the  ICPs  to  compute  the  reorder 
quantity,  Q,  and  the  reorder  point,  R,  for  consumables  will  be  discussed.  The  following 
topics  wll  be  covered: 

1.  Total  variable  cost  equation. 

2.  Model  assumptions. 

3.  Derivations  of  Q  and  R. 

4.  Constraints. 

1.  Total  Variable  Cost  Equation 

The  Department  of  Defense  (DOD)  provides  policies  for  determining  inventory 
levels  in  DOD  Instruction  4140.39  and  states  that  the  objective  is  "to  minimize  the  total 
of  variable  order  and  holding  costs  subject  to  a  constraint  on  time-weighted, 
essentiality-weighted  requisitions  short."  [Ref.  6:  pp.  2] 

The  total  variable  cost  (TVC)  equation  in  DOD  Instruction  4140.39  looks  like: 

TVC  —  Order  Cost  -I-  Inventory  Holding  Cost  +  Shortage  Cost 

The  order  cost,  inventory  holding  cost  and  shortage  cost  are  functions  of  Q  and  R. 

The  UICP  model,  based  on  DOD  Instruction  4140.39  is  the  minimization  of  an  an¬ 
nual  variable  cost  equation  composed  of  the  sum  of  three  terms: 

Ordering  Cost  -  an  average  number  of  orders  per  year  times  the  administrative  cost 
of  placing  an  order  or  contract. 

Holding  Cost  -  average  number  of  units  of  stock  on  hand  at  any  random  point  in 
time,  times  the  cost  to  hold  one  unit  of  stock  in  inventory  for  a  year. 

Shortage  Cost  -  average  number  of  requisitions  on  backorder  at  any  random  point 
in  time,  times  the  cost  of  not  filling  one  requisition  for  a  year  times  the 
militaiy-  essentiality  of  the  item.  [Ref.  6:  pp.  3-A-3] 

The  total  annual  variable  cost  (TVC)  equation  is: 

n  n  s 

/*1  /ml 

where: 

/  =  item  index; 
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n  =  total  number  of  items  in  a  inventor>^; 

B  =  expected  number  of  units  of  item  i  backordered  at  any  random  point  in  time; 

).  =  shortage  cost  of  one  requisition  backordered  for  one  year; 

E  =  militaiy’  essentiality  of  item  i;  and 
S  =  expected  number  of  units  per  requisition  for  item  i. 

From  reference  6  the  formula  for  B,{Q„  R)  is; 

BiiQi,  Rd  =  -  Ri)[F,{x  +  -  fXx;  Hi)]dx- 

where; 

F,(  • )  =  cumulative  probability  distribution  of  leadtime  demand. 

2.  Model  Assumptions 

The  ICPs  use  the  following  key  simplifying  assumptions  in  the  LTCP  consum¬ 
ables  inventor}'  model.  While  these  assumptions  may  not  reflect  the  "real  world"  envi¬ 
ronment  in  which  the  inventor}’  model  operates,  the  assumptions  are  necessary  to 
streamline  and  simplify  the  computation  of  Q  and  R.  The  assumptions  are; 

1.  Steady  state  demand  environment  exists  (the  demand  distribution  does  not  change 
over  time). 

2.  Infinite  procurement  quantities  are  possible. 

3.  Continual  review  of  the  inventory  position  occurs. 

4.  The  quantitv  is  sufficient  to  increase  the  inventory  position  up  to  a  value  of 
Q  +  R. 

5.  Procurement  cost,  A,  is  a  constant,  independent  of  order  quantity,  Q. 

6.  Demands  and  other  model  parameters  associated  with  different  items  are  inde¬ 
pendent. 

3.  Derivations  of  Q  and  R 

a.  Derivation  of  Order  Quantity  ( Q) 

The  derivation  presented  in  this  section  was  taken  from  reference  8.  The 
optimal  order  quantity,  Q,  is  determined  by  solving; 


where; 
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eB{Q,R) 


oQ 


^~dQ~Q  J +  2;  m)  -  F{x-,  ,j.))dx. 


However,  this  is  an  equation  involving  an  implicit  function  of  Q  and  it  is  difficult  to 

solve  explicitly.  It  can  be  solved  iteratively,  but  those  calculations  would  be  tedious  for 

large  inventory  systems  containing  a  large  number  of  items.  For  this  reason,  the  ICP's 
cTVC 

Ignore  this  term  in  — — —  =  0  and  solve  for  Q  using  the  remaining  terms.  The  result  is 


SQ 


Qeoq- 


%DA 
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b.  Derivation  of  Reorder  Point  (R) 

The  optimal  reorder  point,  R,  is  determined  from 


cR 


cR 


Limiting  arguments  to  eliminate  Q  from  this  formula  results  in  an  optimal  R  satisfying 

SIC 


R1SK=  l-F{R\n)  = 


SIC  +  aIVE  ’ 


where: 

W  =  expected  requisitions  per  quarter. 

4.  Constraints 

The  LTCP  consumable  inventory  model  imposes  several  constraicr ;  on  the  op¬ 
timal  Q  and  R  from  the  total  variable  cost  equations.  These  constraints  are  imposed 
because  many  "real  world"  factors  cannot  be  included  in  the  cost  equations.  These 
constraints  will  be  discussed  in  this  section. 
a.  Reorder  Quantity  Constraints 

The  ICP's  impose  constraints  on  the  quantity  they  actually  order  so  that 

A 

their  order  quantity,  Q,  is  determined  from  the  following  formula: 

A 

Q  =  mmmum{\2D,  maximum{Q^QQ,  I,  D)}. 

This  formula  insures  that: 

1 .  The  reorder  quantity  is  at  least  1 . 
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2.  An  item  is  reordered  no  more  than  once  per  quarter  (by  ordering  at  least  D,  the 
quarterly  demand). 

3.  The  maximum  size  of  an  order  is  no  more  than  3  years'  demand  (12D)  per 
DODINST  4140.39. 

b.  Reorder  Point  Constraints 

The  reorder  point  is  constrained  in  two  steps.  First,  the  right  hand  side  of 
the  RISK  formula  is  computed.  That  RISK  value  is  then  constrained  to  a  value  of  no 
less  than  an  ICP-set  minimum  risk  and  no  greater  than  an  ICP-set  maximum  risk.  These 
risk  limits  are  established  to  prevent  under  or  over  investment  in  any  given  item. 

The  reorder  point,  R*,  for  the  constrained  risk  is  computed  next.  The  dis¬ 
tribution  used  for  any  item  is  dependent  upon  the  demand  characteristics  of  that  item. 
LTCP  limits  the  choices  of  distributions  to  the  normal,  Poisson  and  negative  binomial. 
The  choice  is  based  on  a  comparison  of  /i,  where  n  =  D  x  L,  to  an  ICP-set  parameter, 
called  the  probability  .breakpoint.  If  is  less  than  the  probability  breakpoint  a  negative 
binomial  probability  distribution  is  used.  If  ft  is  greater  than  or  equal  to  the  probability 
breakpoint  a  normal  probability  distribution  is  selected.  A  Poisson  probability  distrib¬ 
ution  is  used  for  veiy  low-demand  items,  specially  designated  by  a  "Mark  Code  of  0". 

The  reorder  point,  R*,  is  then  further  constrained  to  no  less  than  zero  or 
an  ICP-set  value,  called  the  Numeric  Stockage  Objective  (NSO  is  a  minimum  stockage 
icvel  established  by  an  ICP  on  low  demand  items  for  insurance  stockage  purposes),  no 
less  than  an  ICP-set  percentage  of  ft  or  the  number  of  wholesale  stockpoints  carrying  the 
wholesale  level  for  the  item  (called  policy  receivers)  designated  by  the  ICP,  and  no  larger 
than  an  ICP-set  percentage  of  ft  with  on-hand  assets  not  exceeding  the  shelf  life  quan- 

A 

tity.  The  result  is  the  final  constrained  reorder  point,  R  . 

5.  Selection  of  2 

The  shortage  parameter,  2,  is  chosen  so  that  the  resulting  inventor^’  levels  meet 
the  S.MA  goal  of  85°/'o.  Later  it  is  reduced  if  sufficient  funds  are  not  made  available  (not 
a  current  problem  under  stock  funding). 
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III.  MEAN  SUPPLY  RESPONSE  TIME  MODEL 

This  chapter  explains  both  the  theory  and  the  computational  procedures  used  in  the 
mean  supply  response  time  (MSRT)  wholesale  consumable  inventory  model. 

A.  THE  MSRT  MODEL 

In  the  MSRT  model  the  objective  is  to  determine  the  depth  of  each  item  which 
minimizes  the  aggregate  mean  supply  response  time  over  all  items  within  a  given  four¬ 
digit  cog,  subject  to  a  budget  constraint  or  maximum  investment  level  for  that  four-digit 
cog.  (As  stated  in  chapter  I,  the  mean  supply  response  time  is  defined  as  the  mean  time 
it  takes  the  supply  system  to  satisfy  the  demand  for  an  item.) 

Mathematically  the  MSRT  model  seeks  the  optimal  maximum  integer  depths, 
SIF,  >  0  (i  =  1,  2, ...,  n),  which  minimize; 


MSRT^ 


n 

X  MSRTiSW^ 

i=\ _ 

n 

/=! 


subject  to: 


n 

^  ^  Investment  Level', 

(=1 


where: 

i  =  item  index; 

D  =  forecasted  average  quarterly  demand  for  item  i; 
SW  =  the  depth  of  item  i;  and 
C  =  unit  price  of  item  i. 


It  has  been  shown  [Ref.  1]  that: 


MSRT, 


TWUS, 

A  ’ 
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where  TWUS  stands  for  the  expected  time-weighted  units  short  per  year.  It  is  defined 
[Ref.  4:  pp.  185]  as: 


SH' 

TPyUS,(S  =  (Mi -s  IV,)  +J^(SfV,-  X,)  X  p,(x,;  m,); 

where  p(  x;  /x  )  is  the  Poisson  or  negative  binomial  probability  mass  function  for  demand 
X  during  lead  time  and  m  =  Dx  L.  If  demand  can  be  approximated  by  the  normal  dis¬ 
tribution  then  the  value  of  p(  x;  )  for  integer  x  values  can  be  obtained  by  using  a 
continuity  correction  on  the  normal.  Therefore,  the  objective  function  can  be  rewritten 
as: 


MSRT^ 


n 

i=\ 


and,  because  the  denominator  is  a  constant,  the  objective  function  reduces  to: 


n 

/=! 


Conceptually,  time- weighted  units  short  for  an  item  are  the  number  of  demands  that 
can't  be  satisfied  by  issues  from  stock,  weighted  by  the  length  of  time  that  each  demand 
remains  unsatisfied  (in  a  backorder  position).  In  the  case  of  consumables  these  demands 
normally  remain  backordered  until  a  procurement  arrives. 

The  constraint  in  the  MSRT  model  has  been  established,  for  comparison  purposes, 
as  the  total  investment  level  resulting  from  the  UlCP  model  for  the  same  four-digit  cog. 

A  A 

That  value  is  obtained  by  adding  the  <2  and  R,  values  to  get  the  maximum  UICP  depth 
for  each  item,  multiplying  this  depth  by  the  unit  cost  C] ,  and  summing  these  products 
over  all  i,  i  =  I,  2, ...,  n. 

B.  MSRT  MODEL  SOLUTION 

The  optimization  technique  used  on  the  MSRT  model  is  marginal  analysis  because 
of  the  integer  nature  of  demand  and  hence  the  depth  SW.  Basic  to  the  marginal  analysis 
approach  is  the  ratio  /?/?„  where: 
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TWUSi{X-l)-TWUS^X) 


This  represents  the  ratio  of  the  improvement  in  TWUSt,  when  the  depth  of  item  i  is  in¬ 
creased  from  X-1  to  X,  to  the  cost  C,  incurred  when  the  depth  is  increased  by  that  one 
unit. 

To  compute  an  optimal  depth  for  the  MSRT  model  using  marginal  analysis,  all 
items  are  initialized  with  a  depth  of  zero  (SW  =  0)  and  RR  is  computed  for  X  =  1  for  all 
i.  That  item  having  the  largest  value  of  RR,-  is  then  increased  by  one  in  its  depth  (SW). 
Its  RR  value  is  next  computed  with  X  =  2  and  all  RR,  values  are  compared  again.  That 
item  whose  RR  is  largest  has  its  SW  increased  by  one  unit,  its  X  value  is  increased  by 
one,  and  a  new  RR  value  is  computed  for  it.  After  each  unit  increase  in  an  item  i,  the 
left-hand  side  of  the  investment  constraint  is  increased  by  C,.  This  procedure  of  com¬ 
paring  RR,  values,  increasing  one  item's  depth  and  reducing  the  budget  constraint  is  re¬ 
peated  until  the  addition  of  the  next  unit  would  cause  the  cumulative  investment  level 
for  the  items  within  that  four-digit  cog  to  exceed  the  investment  level  the  LTCP  model 
produced  for  the  same  four-digit  cog.  At  this  point  the  calculations  are  terminated.  The 
calculations  could  be  continued,  however,  for  only  those  items  whose  unit  prices  remain 
low  enough  not  to  exceed  the  budget  constraint  on  the  next  step. 

Each  time  an  item  goes  through  the  procedure  of  computing  the  improvement 
brought  about  by  adding  one  unit  to  its  depth,  its  mean  supply  response  time  in  days  is 
also  computed  using  the  following  formula: 


MSRTi  = 


91  X  rifX’sxsn’,.) 
Di 


If  the  item's  MSRT  is  less  than  0.01  it  is  excluded  from  further  depth  increases.  This 
value  was  set  as  a  lower  bound  on  any  item's  MSRT  since  further  reductions  would  not 
significantly  enhance  system  performance.  The  value  of  0.01  is,  however,  arbitrary. 

As  mentioned  earlier,  to  produce  an  optimal  depth,  SW,  using  the  MSRT  model  the 
depth  of  each  item  needs  to  be  initialized  at  zero.  If  initialized  at  zero  the  model 
produced  a  depth  of  zero  for  many  low  demand  items  used  in  the  comparative  analyses. 
Although  this  is  correct,  depths  of  zero  are  rarely  allowed  in  the  LTCP  model.  In  the 

A  A 

LTCP  model  the  Q  +  R  depth  is  almost  always  one  or  greater  because  the  order  quantity 
is  almost  always  constrained  to  be  one  and  the  reorder  point  is  constrained  to  no  less 
than  zero.  Despite  the  constraints  in  the  LTCP  model  for  very  slow  movers  this  is  not 
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an  accurate  representation  of  the  way  these  items  are  actually  managed.  Usually  these 
items  have  no  material  on-hand  and  no  material  on-order  until  a  requisition  is  received. 

A  A 

The  true  reorder  point  then  is  -1,  the  reorder  quantity  is  1  and  the  0+  depth  is  zero. 

In  this  research  effort  the  MSRT  model  was  run  under  three  different  conditions: 
depth  initialized  at  zero;  depth  initialized  at  one;  and  depth  initialized  at  the  LTCP 
computed  reorder  point.  With  the  depth  initialized  at  one  the  results  may  not  be  opti¬ 
mal,  but  they  provide  a  basis  for  a  comparison  with  the  existing  LTCP  model. 

The  MSRT  model  w’as  run  with  the  depth  initialized  at  the  reorder  point  computed 
by  the  LTCP  model  to  assess  the  impact  on  the  different  measures  of  effectiveness  and 
computer  run  times  of  never  having  a  depths  less  than  the  LTCP  reorder  points  com¬ 
puted  by  LTCP. 
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IV.  SOFTWARE,  HARDWARE  AND  DATA 

This  chapter  describes  the  data  that  was  used  to  run  the  models  and  the  computer 
hardware  and  software  utilized. 

A.  SOFTWARE  AND  HARDWARE 

The  original  programming  of  the  MSRT  model  was  done  in  FORTRAN  66  by 
Professor  Alan  McMasters.  The  first  step  in  this  thesis  research  was  to  convert  the  or¬ 
iginal  programming  to  WATFOR  77.  The  conversion  was  made  to  allow  for  easier 
program  formulation,  analysis  (debugging)  and  structuring.  The  completed  program  is 
in  WATFOR  77  Release  2.5.  The  operating  system  was  VM  CMS.  The  hardware  on 
which  the  programming  and  running  of  the  models  was  done  was  the  IBM  3033  model 
mainframe  at  the  W.R.  Church  Computer  Center  at  the  Naval  Postgraduate  School, 
Monterey,  California. 

B.  DATA 

The  data  that  was  used  for  the  comparative  analyses  was  from  a  tape  containing  the 
data  for  all  the  IH  cog  National  Item  Identification  Numbers  (NTINs)  on  file  at  SPCC 
in  the  spring  of  19S5.  The  tape  contained  113,647  IH  records  with  each  record  pro¬ 
viding  information  on  one  NTIN.  The  records  were  in  the  CARES  (Computation  and 
Research  Evaluation  System)  format. 

Because  of  the  program  size  and  the  limitation  on  virtual  storage  and  CPC  time 
available  for  this  research,  only  select  data  w'as  used.  Computer  runs  were  restricted  to 
one  hour  of  CPU  time  due  to  policies  of  the  NPS  Computer  Center.  Data  samples  were 
therefore  selectively  restricted  in  size  to  enable  program  processing  within  the  one  hour 
time  limit.  Data  samples  were  also  restricted  in  size  because  of  the  large  number  of  ar¬ 
rays  used  in  the  program  and  the  limited  virtual  storage  available  to  an  individual  user. 

The  process  used  for  selecting  data  was  to  first  compute  the  average  quarterly  de¬ 
mand  for  all  items  within  each  four-digit  cog.  Next,  enough  four-digit  cogs  were  selected 
to  represent  the  entire  range  of  average  quarterly  demands.  Finally,  if  all  the  NTINs 
within  a  given  four-digit  cog  couldn't  be  processed  through  the  models  within  time  and 
space  limitations,  a  subset  of  those  items  was  selected  which  had  an  average  demand 
closely  approximating  that  of  the  four-digit  cog's  entire  population. 
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Table  1  is  a  breakdown  of  the  entire  IH  data  set  and  information  on  the  data  sets 
that  were  used  for  the  comparisons.  The  first  entry  in  Table  1  shows  the  items  from  the 
data  that  were  not  assigned  a  four-digit  cog. 


Table  I.  BREAKDOWN  BY  4-DIGIT  COG  OF  THE  DATA  BASE 


COG 


IH 


IHM 


1HN2 


1HN3 


1HN4 


IHSl 


1HS2 


IHS3 


1HS4 


IHOA 


IHOB 


IHOC 


IHOD 


IHOE 


IHOF 


IHIA 


IHIB 


IHIC 


IHID 


IHIE 


IHIF 


1H2A 


1H2B 


1H2C 


1H2D 


1H2E 


1H2F 


1H.3A 


1H3B 


1H3C 


1H3D 


1H3E 


1H3F 


1H4A 


1H4B 


1H4C 


1H4D 


1H4E 


1H4F 


TOTAL 
NUMBER  OF 
ITE.MS 


1566 


1026 


1893 


13941 


568 


983 


11700 


11 


13 


144 


818 


2005 


18900 


71 


174 


2037 


582 


1091 


4711 


79 


32 


613 


1851 


3150 


18952 


AVER.AGE 

QUARTERLY 

DEMAND 


0.81 


325.52 


25.89 


1.25 


0.40 


86.16 


21.90 


1.13 


2.43 


526.71 


22.37 


74.42 


4567.95 


3.50 


0.88 


232.50 


12.87 


0.94 


41.13 


11.18 


0.39 


157.66 


25.00 


1.43 


236.03 


4.02 


71.16 


11.08 


1.28 


78.55 


8.57 


1.69 


98.23 


5.73 


0.62 


167.26 


5.61 


0.84 
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V.  MEASURES  OF  EFFECTIVENESS 


This  chapter  discusses  the  Measures  of  Effectiveness  (MOEs)  used  in  the  compar¬ 
ative  evaluation  and  gives  the  computer  resul*s  for  each.  The  following  is  a  list  of  MOEs 
used: 

1.  Mean  Supply  Response  Time  (MSRT). 

2.  Supply  Material  Availability  (SMA). 

3.  Average  Days  Delay  (ADD). 

4.  Average  Days  Delay  for  Backordered  Requisitions  (ADDBO). 

5.  Days  Of  Safety  Level  (DOST). 

6.  CPU  Run  Time. 

7.  Depth  Churn. 

SMA  (specifically  SMA2),  ADD  and  ADDBO  are  standard  CARES  (Computation  and 
Research  Evaluation  System)  .MOEs  [Ref.  9]. 

A.  MEAN  SUPPLY  RESPONSE  TIME  (MSRT) 

Mean  Supply  Response  Time  is  the  mean  time  in  days  taken  by  the  supply  system 
to  respond  to  the  demand  for  an  item.  The  aggregate  MSRT  for  each  four-digit  cog  is 
computed  as  follows: 


MSRT= 


^[91  X  riTUS,] 

f=l _ 

n 


1=1 


where: 

TWU  S  =  Time-weighted  units  short  in  days  for  item  i; 

D  =  Forecasted  average  quarterly  demand  for  item  i;  and 
n  =  Number  of  items  in  the  four-digit  cog. 

Table  2  gives  the  results  for  the  MSRT  measurements  for  the  MSRT  and  LTCP  models 
with  the  MSRT  model  depth  initialized  (ID)  at  zero,  one  and  the  LTCP  reorder  point. 


19 


Table  2.  MODEL  COMPARISONS  FOR  MSRT  (IN  DAYS) 


COG 

SAMPLE 

AVERAGE 

QUARTERLY 

DEMAND 

UICP 

MODEL 

MSRT 

MODEL 

ID  =  0 

MSRT 

MODEL 

ID  =  1 

MSRT 

MODEL 

ID  =  .R 

IHOD 

4567.95 

1.384 

0.394 

0.401 

0.606 

IHOA 

511.46 

1.659 

0.369 

0.369 

0.403 

1H2D 

236.03 

1.189 

0.653 

0.653 

0.671 

1H2A 

140.91 

1.920 

1.221 

1.221 

1.242 

IHSl 

86.16 

1.096 

0.491 

0.491 

0.531 

1H2B 

24.17 

4.192 

1.824 

1.852 

2.217 

1H3E 

8.57 

12.958 

4.209 

4.219 

4.749 

1H4F 

0.84 

36.871 

3.412 

5.510 

6.877 

IHIF 

0.40 

24.738 

0.027 

0.455 

0.496 

1.  MSRT  Analysis  and  Observations 

As  would  be  expected,  given  its  objective  function  of  minimizing  mean  supply 
response  time,  the  MSRT  model  provides  better  mean  supply  response  time  than  the 
LTCP  model  in  all  cases.  The  improvement  in  performance  becomes  greater  as  the  av¬ 
erage  demand  of  the  population  decreases. 

The  MSRT  model  with  depth  initialized  at  zero  provides  the  best  results,  but  the 
difference  in  performance  between  the  three  forms  of  the  MSRT  model  is  negligible  until 
the  average  quarterly  demand  of  the  sample  gets  very  small. 

B.  SUPPLY  MATERIAL  AVAILABILITY  (SMA) 

Supply  Material  Availability  is  the  percent  of  all  requisitions  that  are  satisfied  from 
on-hand  stock  over  a  year.  The  Chief  of  Naval  Operations  has  set  the  current  SMA  goal 
at  85%.  For  purposes  of  this  research  SMA  is  computed  two  different  ways,  called 
SMAl  and  SMA2.  SMAl  is  a  demand-based  measure  calculated  as  one  minus  the  ex¬ 
pected  number  of  backorders  divided  by  the  average  quarterly  demand.  SMA2  is  a 
requisition-based  measure  calculated  as  one  minus  the  number  of  requisitions  short  per 
year  divided  by  the  average  annual  requisition  frequency.  The  formulas  for  each  are  as 
follows: 
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=  1  - 


t=l _ 

n 

/=1 


SMA2  =  1  - 


n 

^[/?spy,] 

/=! _ 

n 

AxY^m 

1=1 


where: 

EBO  =  Expected  number  of  backorders  per  year; 

RSPY  =  Requisitions  short  per  year;  and 
W  =  Quarterly  requisition  frequency. 

Tables  3,  4,  and  5  present  the  results  for  these  two  SMA  MOEs  for  the  MSRT 
model,  with  the  depth  for  each  item  initialized  at  0,  1  and  the  reorder  point  (RP),  re¬ 
spectively. 


Table  3.  MODEL  COMPARISON  FOR  SMA  WITH  THE  MSRT  MODEL 


HAVING  ID  =  0 


COG 

SAMPLE 

AVG. 

DEMAND 

UICP 

SMAl 

MSRT 

SMAl 

UICP 

SMA2 

MSRT 

SMA2 

IHOD 

4567.95 

99.98 

99.98 

92.75 

87.38 

IHOA 

511.46 

99.86 

99.87 

97.67 

96.92 

1H2D 

236.03 

99.73 

99.79 

94.60 

95.51 

1H2A 

140.91 

99.54 

99.61 

94.40 

94.86 

IHSl 

86.16 

99.47 

99.69 

96.40 

97.40 

1H2B 

24.17 

98.62 

99.32 

88.90 

92.75 

1H3E 

8.57 

97.54 

98.62 

91.95 

94.37 

1H4F 

0.84 

93.30 

99.31 

86.56 

98.74 

IHIF 

0.40 

84.35 

99.66 

84.36 

99.03 
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Table  4.  MODEL  COMPARISON  FOR  SMA  WITH  THE  MSRT  MODEL 
HAVING  ID  =  1 


COG 

SAMPLE 

AVG. 

DEMAND 

UICP 

SMAl 

MSRT 

SMAl 

UICP 

SMA2 

MSRT 

SMA2 

IHOD 

4567.95 

99.98 

99.98 

92.75 

87.32 

IHOA 

511.46 

99.86 

99.87 

97.67 

96.92 

1H2D 

236.03 

99.73 

99.79 

94.60 

95.51 

1H2A 

140.91 

99.54 

99.61 

94.40 

94.86 

IHSl 

86.16 

99.47 

99.69 

96.40 

97.39 

1H2B 

24.17 

98.62 

99.31 

88.90 

92.62 

1H3E 

8.57 

97.54 

98.61 

91.95 

94.33 

1H4F 

0.84 

93.30 

99.00 

86.56 

98.31 

IHIF 

0.40 

84.35 

99.02 

84.36 

98.75 

Table  5.  MODEL  COMPARISONS  FOR  SMA  WITH  THE  MSRT 
MODEL  HAVING  ID  =  RP 


COG 

SAMPLE 

AVG. 

DEMAND 

UICP 

SMAl 

MSRT 

SMAl 

UICP 

SMA2 

MSRT 

SMA2 

IHOD 

4567.95 

99.98 

99.98 

92.75 

89.09 

IHOA 

511.46 

99.86 

99.86 

97.67 

96.73 

1H2D 

236.03 

99.73 

99.78 

94.60 

95.43 

IH2A 

140.91 

99.54 

99.60 

94.40 

94.55 

IHSl 

86.16 

99.47 

99.64 

96.40 

96.97 

1H2B 

24.17 

98.62 

99.16 

88.90 

91.11 

1H3E 

8.57 

97.54 

98.32 

91.95 

93.26 

1H4F 

0.84 

93.30 

98.46 

86.56 

97.42 

IHIF 

0.40 

84.35 

98.64 

84.36 

98.66 

1.  SMA  Analysis  and  Observations 

Relative  to  the  LTCP  model,  the  MSRT  model  provides  improved  or  equivalent 
SMA  performance  in  all  cases  except  for  SMA2  in  the  case  of  very  high  demand  items. 
For  both  measures  of  SMA  the  improvement  in  performance  provided  by  the  MSRT 
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model  becomes  greater  as  the  average  demand  of  the  group  of  items  processed  becomes 
smaller.  It  is  interesting  to  note  that,  in  general,  the  performance  of  both  the  UlCP  and 
MSRT  models  deteriorates  as  the  average  demand  of  the  population  decreases.  The 
better  performance  of  the  MSRT  model  then  results  from  the  fact  that  its  SMA  values 
decrease  at  a  slower  rate  than  those  of  the  LTCP  model. 

Comparison  of  the  performance  of  the  three  MSRT  models,  (depth  initialized 
at  0,  I  or  the  reorder  point)  shows  remarkably  similar  performance  with  very  few  differ¬ 
ences  in  SMA  greater  than  ±  1  %.  Again,  the  MSRT  model  with  depth  initialized  at 
zero  provides  the  best  results. 

C.  AVERAGE  DAYS  DELAY  (ADD) 

Average  Days  Delay  measures  how  long  on  the  average  the  supply  system  takes  to 
fill  a  requisition,  on  the  average  aggregated  over  all  requisitions.  ADD  for  each  four¬ 
digit  cog  is  computed  two  different  ways  as  follows: 


ADD\  = 


ri 

365  X  YjTWEB^ 


/=! 


n 

4x^[in 

/asl 


ADD2 


w'here: 

TWEB  =  Time-weighted  average  number  of  requisitions  backordered  per  year  for 
item  i. 

ADDl  is  strictly  a  requisition-based  calculation.  ADD2  begins  as  requisition-based 
but  is  converted  to  demand-based  by  factoring  in  average  requisition  size.  Tables  6  and 
7  give  the  ADD  results. 
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Table  6.  AVERAGE  DAYS  DELAY  (ADDl) 


COG 

SAMPLE 

AVERAGE 

DEMAND 

UlCP 

ADDl 

MSRT 

ADDl 

ID  =  0 

MSRT 

ADDl 

ID  =  1 

MSRT 

ADDl 

ID  =  ^ 

I  HOD 

4567.95 

3.80 

27.69 

27.79 

16.11 

I  HO  A 

511.46 

1.86 

4.80 

4.80 

3.72 

1H2D 

236.03 

7.92 

9.50 

9.50 

8.63 

1H2A 

140.91 

7.82 

8.70 

8.70 

8.94 

IHSl 

86.16 

3.77 

3.95 

3.98 

4.16 

1H2B 

24.17 

26.94 

18.99 

19.44 

23.62 

1H3E 

8.57 

21.39 

17.81 

18.00 

20.27 

1H4F 

0.84 

65.15 

4.79 

6.85 

11.49 

IHIF 

0.40 

122.74 

7.95 

10.58 

11.35 

Table  7.  AVERAGE  DAYS  DELAY  (ADD2) 


COG 

SAMPLE 

AVER.^GE 

DEMAND 

UICP 

ADD2 

MSRT 

ADD2 

ID  =  0 

MSRT 

ADD2 

ID  =  1 

MSRT 

ADD2 

ID  =  R 

IHOD 

4567.95 

2.03 

0.58 

0.59 

0.89 

IHOA 

511.46 

2.27 

0.64 

0.64 

0.71 

1H2D 

236.03 

1.94 

1.03 

1.03 

1.05 

1H2A 

140.91 

2.52 

1.45 

1.45 

1.49 

IHSl 

86.16 

2.21 

1.28 

1.28 

1.37 

1H2B 

24.17 

13.10 

9.22 

9.30 

10.25 

1H3E 

8.57 

16.07 

5.30 

5.35 

6.19 

1H4F 

0.84 

50.44 

6.31 

8.55 

11.26 

IHIF 

0.40 

82.66 

26.46 

33.18 

35.22 

1.  ADD  Analysis  and  Observations 

The  results  of  the  ADDl  MOE  are  slightly  inconsistent.  In  the  case  of  ADDl 
the  MSRT  model  performs  worse  than  the  LTCP  model  for  high  demand  items  and 
better  than  the  LTCP  model  for  low’  demand  items.  Without  exception,  however,  the 
.MSRT  models  perform  better  than  the  LTCP  model  for  ADD2.  Once  again,  the  im¬ 
provement  in  performance  increases  as  the  average  demand  of  the  group  decreases.  The 


24 


performances  of  the  three  forms  of  the  MSRT  models  are  similar  with  the  depth  initial¬ 
ized  at  zero  being  the  best. 

D.  AVERAGE  DAYS  DELAY  FOR  BACKORDERED  REQUISITIONS  (ADDBO) 
Average  Days  Delay  for  Backordered  Requisitions  measures  how  long  an  average 
backordered  requisition  remains  on  backorder  before  it  is  filled.  ADDBO  is  a  standard 
CARES  requisition-based  .MOE.  The  formula  is: 

Table  8  gives  the  ADDBO  results. 


Table  8.  AVERAGE  DAYS  DELAY  FOR  BACKORDERED  REQUISITIONS 


(ADDBO) 


COG 

SAMPLE 

AVERAGE 

DEMAND 

UICP 

ADDBO 

MSRT 

ADDBO 

ID  =  0 

MSRT 

ADDBO 

ID  =  1 

MSRT 

ADDBO 

ID  =  R 

IHOD 

4567.95 

52.35 

219.43 

219.11 

147.63 

IHOA 

511.46 

80.04 

155.82 

155.82 

113.86 

IH2D 

236.03 

146.64 

211.63 

211.63 

188.77 

1H2A 

140.91 

139.60 

169.08 

169.08 

163.99 

IHSl 

86.16 

104.61 

152.12 

152.28 

137.50 

1H2B 

24.17 

242.62 

262.20 

263.23 

265.82 

1H3E 

8.57 

265.76 

316.39 

316.95 

300.76 

1H4F 

0.84 

484.59 

379.26 

405.95 

446.11 

IHIF 

0.40 

784.95 

822.79 

845.48 

845.73 

1.  ADDBO  Analysis  and  Observations 

An  quick  analj  sis  of  the  results  for  this  MOE  demonstrates  a  clear  advantage 
to  the  UICP  model. 

E.  DAYS  OF  SAFETY  LEVEL 

Days  of  safety  level  measures  the  amount  of  average  stock  on  hand  to  protect 
against  a  stockout  during  leadtime.  The  higher  the  days  of  safety  level  the  lower  the  risk 
of  a  stockout.  The  formula  for  days  of  safety  level,  which  is  computed  by  four-digit  cog, 
is: 
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365  X  Yj[Q  X  5L,1 

DOSL^—T^ - ; 

£[4xqxZ),] 

i=l 

where: 

DOSL  =  Cost  weighted  average  safety  level  in  days  of  supply  at  the  forecasted  mean 
quarterly  demand  rate  of  usage; 

SL  =  Safety  Level  (Item  maximum  depth  —  reorder  quantity  —  leadtime  demandin 
units  for  item  i);  and 

C  =  Unit  price  for  item  i. 

THe  UICP  reorder  quantity  was  also  used  for  the  reorder  quantity  of  the  MSRT  model. 
Table  9  gives  the  days  of  safety  level  results. 


Table  9.  DAYS  OF  SAFETY  LEVEL  (DOSL) 


COG 

SAMPLE 

AVERAGE 

DEMAND 

UICP 

DOSL 

MSRT 

DOSL 

ID  =  0 

MSRT 

DOSL 

ID  =  1 

MSRT 

DOSL 

ID  =  R 

IHOD 

4567.95 

48.03 

69.38 

68.65 

51.81 

IHOA 

511.46 

138.07 

146.81 

146.81 

138.37 

1H2D 

236.03 

222.76 

232.77 

232.77 

222.72 

1H2A 

140.91 

222.25 

224.41 

224.41 

222.24 

IHSl 

86.16 

223.56 

247.14 

246.37 

224.34 

1H2B 

24.17 

164.56 

243.19 

238.35 

190.46 

1H3E 

8.57 

287.67 

360.96 

359.49 

302.76 

1H4F 

0.84 

360.06 

779.76 

649.53 

489.00 

IHIF 

0.40 

413.09 

1181.84 

627.79 

557.59 

1.  Days  of  Safety  Level  Analysis  and  Observations 

An  analysis  of  the  days  of  safety  level  MOE  show’s  that  the  MSRT  model  with 
item  denths  initialized  at  either  0,  1  or  the  reorder  point  clearly  outperforms  the  UICP 
model.  As  before,  the  degree  of  improvement  increases  as  the  average  demand  of  the 
group  decreases. 
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F.  CPU  RUN  TIME 


Central  Processor  Unit  (CPU)  run  time  measures  how  long  the  computer  takes  to 
compute  the  optimal  depth  for  each  model.  CPU  run  time  is  an  important  MOE  be¬ 
cause  any  additional  run  time  may  affect  the  ICP's  if  their  ADP  capacity  is  already 
nearing  its  limits.  Table  10  gives  the  CPU  run  time  results. 


Table  10.  CPU  RUN  TIME  (CPURJ)  IN  MINUTES 


COG 

SA.MPLE 

AVERAGE 

DEMAND 

SAMPLE 

SIZE 

UICP 

CPURT 

MSRT 
CPURT 
ID  =  0 

.MSRT 
CPURT 
ID  =  1 

MSRT 

CPURT 

ID=  R 

IHOD 

4567.95 

11 

0.01 

14.18 

13.84 

5.92 

IHOA 

511.46 

100 

0.10 

38.60 

37.96 

16.69 

1H2D 

236.03 

79 

0.01 

14.22 

14.07 

6.58 

1H2A 

140.91 

100 

0.02 

15.82 

15.63 

5.99 

IHSl 

86.16 

282 

0.05 

49.13 

47.85 

29.00 

1H2B 

24.17 

492 

0.08 

48.71 

45.51 

26.73 

1H3E 

8.57 

583 

0.09 

29.80 

28.48 

13.30 

1H4F 

0.84 

882 

0.15 

9.05 

8.06 

5.06 

IHIF 

0.40 

1420 

0.20 

8.49 

8.05 

5.25 

1.  CPURT  Analysis  and  Obsenations 

As  expected,  the  MSRT  model  with  its  marginal  analysis  procedure  for  deter¬ 
mining  optimal  item  depths  under  a  budget  constraint  requires  significantly  longer  CPU 
run  time  than  the  LTCP  model.  Initializing  the  depth  of  each  item  at  the  UICP  com¬ 
puted  reorder  point  decreased  the  MSRT  model  run  time  by  about  50%. 

The  run  times  encountered  for  the  MSRT  model  during  this  research  were  sig¬ 
nificant.  However,  these  can  be  reduced  by  use  of  VS  FORTRAN,  more  efficient  pro¬ 
gramming  techniques  and  the  more  up-to-date  hardware  available  at  the  I  CPs.  The  run 
times  are  probably  comparable  to  those  encountered  when  a  new  A  value  is  being 
computed  selected  by  CARES  for  the  UICP  model.  The  value  of  A  is  not  changed  very 
often.  The  long  runs  to  recompute  maximum  depths  for  the  MSRT  model  need  not 
occur  very  often  either.  Once  the  depths  have  been  computed,  each  item's  MSRT  value 
could  be  used  and  is  available  as  a  subsequent  goal  for  fine  tuning  its  depth  on  a  quar¬ 
terly  basis,  if  that  is  desired. 
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G.  DEPTH  CHURN 


To  assess  the  changes  in  maximum  depth  that  would  be  expected  from  implementing 
the  model,  the  "depth  churn"  between  the  UICP  model  and  the  MSRT  model  with  ID 
=  1  was  computed.  For  item  i,  this  chum  is  defined  as: 

Churrii  —  UICP  Model  Maximum  Depthi  —  MSRT  Model  Depthi 


The  chum  data  in  Appendix  I  represents  the  four-digit  cogs  with  the  most  extreme  and 
least  extreme  churn.  Each  churn  value  is  presented  with  its  frequency  of  occurrence. 
Cogs  IHOA  and  IHOD  showed  the  most  extreme  chum  while  cogs  IH3E  and  1H4F 
showed  the  least. 

1.  Depth  Chum  Analysis  and  Observations 

The  characteristics  of  the  items  that  showed  the  greatest  degree  of  churn,  one 
+  and  one  — ,  for  each  of  the  four-digit  cogs  listed  in  Appendix  1  are  presented  in  Tables 

A 

1 1  and  12  below.  The  Q  value,  which  is  the  same  for  both  models,  is  included  along  with 
each  item's  unit  price  and  quarterly  demand.  Depth  is  in  units  of  inventory. 


Table  II.  COGS  WITH  LEAST  CHURN 


COG 

UICP 

DEPTH 

.MSRT 

DEPTH 

CHURN 

Q 

UNIT 

PRICE 

QTRLY 

DMD 

1H3E 

8151 

10045 

-1894 

805 

38.00 

508.64 

IH3E 

87 

67 

20 

20 

5920.00 

6.19 

1H4F 

984 

1499 

-515 

676 

6.30 

56.24 

1H4F 

40 

25 

15 

20 

7420.00 

1.61 

Keeping  in  mind  that  a  negative  chum  indicates  a  greater  depth  for  the  MSRT 
model,  it  is  clearly  seen  that  the  MSRT  model  places  greater  emphasis  on  low-cost, 
high-demand  items.  This  is  the  expected  result  given  that  the  MSRT  model  criteria  for 
increasing  an  item's  depth  is  based  on  a  comparison  of  the  cost-weighted  ratios  of  im¬ 
provement  in  time-weighted  units  shon. 
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Table  12.  COGS  WITH  GREATEST  CHURN 


COG 

UlCP 

DEPTH 

MSRT 

DEPTH 

CHURN 

A 

Q 

UNIT 

PRICE 

QTRLY 

DMD 

IHOA 

65653 

76059 

-10406 

10411 

4.70 

10078.11 

IHOA 

201 S 

1886 

132 

776 

25.50 

349.17 

IHOD 

50199 

57416 

-7217 

11219 

1.50 

7666.90 

I  HOD 

1294 

1148 

146 

589 

43.50 

288.37 
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VI.  SUMMARY  AND  RECOMMENDATIONS 

A.  SUMMARY 

This  thesis  is  part  of  the  Naval  Postgraduate  School's  continuing  effort  to  introduce 
more  effective  inventory  models  into  the  Naval  Supply  System's  management  of  whole¬ 
sale  level  inventories.  The  specific  objective  of  this  thesis  was  to  develop  an  inventor}^ 
mod’l  for  consumable  items  having  the  objective  function  of  minimizing  the  aggregate 
mean  supply  response  time  (MSRT)  and  then  to  compare  the  results  of  this  model  with 
the  current  consumable  LTCP  model.  The  MSRT  consumable  model  has  the  same 
readiness  related  objective  as  the  recently  implemented  wholesale  provisioning  model  for 
consumables. 

This  objective  was  accomplished  by  developing  the  mathematics  for  the  new  model, 
developing  computer  programs  for  both  models,  processing  ICP  inventory’  management 
data  through  both  models  and  conducting  a  comparative  evaluation  of  the  results  using 
various  measures  of  effectiveness.  The  main  FORTRAN  computer  program  is  presented 
in  Appendix  A.  The  major  subroutines  are  presented  in  Appendices  B  through  H. 

B.  CONCLUSIONS 

This  research  demonstrated  conclusively  that  the  .MSRT  model  improves  supply 
system  performance  for  consumable  items  with  no  increase  in  investment. 

1.  Standard  Measures  of  Effectiveness  Summary 
In  addition  to  being  superior  relative  to  the  aggregate  .MSRT  measure  of  effectiveness, 
the  .MSRT  model  provides  significant  improvements  over  the  LTCP  model  in  the  fol¬ 
lowing  supply  system  performance  measures  for  items  with  a  low,  average  quarterly  de¬ 
mand: 

1.  Supply  material  availability  (S.MA); 

2.  Average  days  delay  (ADD);  and 

3.  Days  of  safety  level. 

The  improvement  provided  by  the  MSRT  model  decreased  as  the  average  quarterly  de¬ 
mand  of  the  groups  of  items  processed  increased.  Performance  of  the  two  models  was 
equivalent  for  the  items  with  a  high  average  quarterly  demand.  This  trend  was  consist¬ 
ently  observed  across  all  measures  of  effectiveness  with  the  exception  of  Average  Days 
Delay  for  Backorders. 
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2.  CPU  Run  Time  Summary’ 

The  CPU  run  time  required  by  the  MSRT  model  with  its  marginal  analysis  op¬ 
timization  procedure  was  significantly  larger  than  the  run  times  required  by  the  UlCP 
model.  The  amount  of  run  time  required  for  any  group  of  items  was  found  to  be  sensi¬ 
tive  to  the  average  quarterly  demand  of  the  four-digit  cog's  population  and  the  number 
of  items  in  the  cog.  Since  the  focus  of  this  research  was  an  initial  comparison  of  the 
performance  of  the  two  inventory  models,  programming  efficiencies  were  not  pursued. 

Reductions  in  CPU  run  time  could  be  obtained  by: 

1.  Using  VS  FORTRAN  instead  of  WATFOR  77. 

2.  Using  more  modem  computer  hardware. 

3.  Using  more  efficient  programming  techniques. 

C.  RECOMMENDATIONS 

Replacement  of  the  UICP  consumable  model  with  the  MSRT  model  should  be 
considered  by  the  ICPs.  While  the  long  computer  run  times  are  a  disadvantage  of  the 
MSRT  model,  it  should  not  be  necessary  to  make  such  runs  every  quarter.  Instead,  each 
item's  depth  can  be  adjusted  to  retain  its  same  MSRT  value  as  when  the  full  four-digit 
cog  optimization  was  last  run.  This  approach  is  similar  to  that  used  to  determine  the  A 
shortage  cost. 

Viewing  consumables  as  a  limiting  case  of  repairables,  research  should  continue  to¬ 
wards  developing  an  .MSRT  model  applicable  to  the  wholesale  inventory  management 
of  repairable  items. 
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APPENDIX  A.  MAIN  PROGRAM 


PROGRAM  REPMOD 

TEST  WHOLESALE  CONSUMABLE  REPLENISHMENT  MODELS 
THIS  PROGRAM  RUNS  BOTH  THE  UICP  MODEL  AND  THE  MSRT  MODEL. 

THE  UICP  MODEL  INITIALIZES  X  AT  ZERO,  ONE  OR  THE  UICP  COMPUTED 
REORDER  POINT  DEPENDING  ON  WHICH  STEPS  ARE  COMMENTED  OUT  IN  THE 
SUBPROGRAM  "MODOPT"  (STEPS  658  TO  670).  THE  PROGRAM  ALSO  COMPUTES 
THE  CHURN  BETWEEN  THE  TWO  MODELS  AND  STORES  IT  IN  A  FILE  THAT  IS 
OPENED  IN  STEP  65. 

CHARACTER*23  NAME2/'MIN  INVESTMENT  LEVEL’/ 

CHARACTER*23  NAME3/ 'MINIMUM  MSRT  ’/ 

CHARACTER*23  NAMEl/'UICP  CONSUMABLES  MODEL  '/ 

CHARACTER*4  Ql/'l  '/,Q2/'4D  '/,Q3/’EOQ  ’/,Q4/'D 

CHARACTER*4  RQl/'l  '/,RQ2/’4D  V,RQ3/’EOQ  '/,RQ4/'D  '/ 

CHARACTER*4  RQ5/ ' UICP' / ,Q5/ ’UICP’ / 

CHARACTER*4  QM,  QMR 
CHARACTER*!  PRIND(1500) 

REAL  BB(10),QBASIC,QTEMP 

INTEGER  STOP( 1500) ,MARK( 1500) ,LOT( 1500) ,MD,PBP( 1500) ,RLC( 1500) 
INTEGER  N,NN,X( 1500) ,NSO( 1500) ,QMIN( 1500) ,NRPR( 1500) ,XMSRT( 1500) 
INTEGER  Y( 1500 ) , QS , Q2 1 , Q2C , Q22 , NPO , NOPT , KK , NI , ROP( 1500 ) 

INTEGER  RESV( 1500) ,BO( 1500) ,PPRLT( 1500) ,PPR1( 1500) ,PPR2( 1500) 
INTEGER  PPR3( 1500) ,PPR4( 1500) ,AD( 1500) ,CD( 1500) ,OHRFI( 1500) 

INTEGER  XCHURN(1500) 

INTEGER  OHNRFK 1500) ,PPR1YR( 1500) ,UPLIM,R0P1,XUICP( 1500) 

REAL  D( 1500) ,PCLT( 1500) ,LAM( 1500) ,SL( 1500) 

REAL  H,RISK,Z( 1500) ,ZN( 1500) ,T( 1500) ,Q( 1500) ,QQ( 1500) 

REAL  Cl( 1500) ,HC( 1500) ,POC( 1500) ,CHRNCST( 1500) ,CHRNTOT 
CHARACTER*4  C0G1( 1500) ,COG2( 1500) 

REAL  MODMST , MODNSF , MSRTG( 10 ) ,MSRTGG , BGT( 10 ) , PVAR( 1500 ) , RST( 1500 ) 
REAL  DMAD( 1500) ,RF( 1500) ,C11( 1500) ,E( 1500) ,RMIN( 1500) ,RMAX( 1500) 
REAL  T0V(2) ,QRI( 1500) ,QRR( 1500) ,QR( 1500) 

REAL  QQR( 1500) ,AS1( 1500) ,SLC( 1500) ,OBSO( 1500) 

CHARACTER*4  C0GG1( 1500) ,COGG2( 1500) 

REAL*8  B.BUDGG 

REAL  RESLT,B19PT1,B19PT2,B19PT3,B19PT4,B19MIN,B19MAX 
COMMON  SN( 1500,9) 

EXTERNAL  MODMST 

C  ***  THE  NEXT  PARAMETER  MUST  BE  SPECIFIED  WHENEVER  A  NEW  COG  IS 
C  ***  INTRODUCED. 

C  ***  THIS  NUMBER  IS  PROVIDED  BY  THE  PRINTOUT  FROM  TEMPDATA  PROGRAM. 
N=1500 

C  ***  THE  NEXT  NUMBER  SPECIFIES  FULL  TABLE  LISTING(NPO*0)  OR  ONLY 
C  ***  A  SUMMARY  (NPO=l) 

NPO=l 

NNN=0 

NN1=0 

C  ***  THE  VALUE  OF  Q  MUST  BE  COMPUTED.  THE  PARAMETER  MQ  TELLS  THE 
C  ***  SUBROUTINE  ORDQAN  WHICH  VALUE  IS  DESIRED;  MQ=1  FOR  Q=l,  MQ=2  FOR 
C  ***  Q=4*D,  MQ=3  FOR  ()=EOQ,  MQ=4  FOR  Q=D,  MQ=5  FOR  QUICP. 

C  ***  QM  IS  EITHER  Ql,  Q2,  Q3,  Q4  OR  Q5,  RESPECTIVELY. 

MQ=5 
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QMR=RQ5 

C  ***  NN  IS  THE  NUMBER  OF  MEASURES  OF 

C  ***  PERFORMANCE;  IF  SMA  AND  MSRT  ARE  USED  THEN  NN=2. 

NOPT=2 

NN=2 

C  ***  A  COUNT  OF  VIABLE  ITEMS  MUST  MADE  (NOT  HAVING  STRANGE  DATA). 

KK=0 

NI=0 

OPEN  (UNIT=10,  FILE='C0A'  STATUS='OLD') 

OPEN(UNIT=11,FILE='CHRNOA' ,STATUS=’NEW' ) 

DO  23  1=1, N 

22  READ(10,898,END=24)COG1(I),COG2(I),MARK(I),(SN(I,J),J=1,9),SLC(I) 
*PRIND(I),OBSO(I),E(I),PCLT(I),NRPR(I),Cl(I),Cll(I),D(I), 
*PVAR(I),RF(I),RESV(I),BO(I),PPRLT(I),PPRl(I),PPR2(I),PPR3(I), 
*PPR4( I) ,OHRFI( I) .OHNRFK I) ,AD( I) ,CD( I) , 
*ASl(I),DMAD(I),LOT(I),NSO(I),QMIN(I) 

898  FORMAT  ( 2A2 , II ,9A1 , lX.Al, 1X,A1 ,6X,F3. 2 ,F3.  3,F4.  2,7X,I4,3F10.  2,10X, 
*  2F10.  2, 218, 515, 17, IX, 17, IX, 218, F8.  0,38X,F10.  2,218,15) 

IF(D(I).GT.  250000.  )GO  TO  22 
IF(D(I).  LT.  0.  0)GO  TO  22 
IF(C1( I).  GT.  999999.  )G0  TO  22 
IF  (C1(I).  LT.  .  01)GO  TO  22 
IF(LOT(I).NE.  0)G0  TO  22 
IF(C0G1(I).NE. '1H')G0  TO  22 
KK=KK+1 
NI=NI+1 

23  CONTINUE 

24  DO  5  1=1, NI 
C0GG1(I)=C0G1(I) 

C0GG2(I)=C0G2(I) 

RMIN(I)=. 10 

IF(COG2(I).EQ.  ’4A'.0R.  COG2(I).EQ.  '3A’.OR.  C0G2(I).EQ.  '2A'.0R. 
*COG2(I).EQ.  '1A'.0R.C0G2(I).EQ.  '0A’)THEN 
RMAX(I)=.  30 
LAM(I)=1500 
PBP(I)=0 
RLC(I)=1 

ELSEIFC  COG2( I ) .  EQ.  ' 4B ’ .  OR.  COG2( I ) .  EQ.  ' 3B ’ .  OR.  COG2( I ) .  EQ.  ' 2B ' .  OR. 
*C0G2(I).EQ.  'IB’.  OR.  C0G2(I).EQ.  'OB')THEN 
RMAX(I)=.  40 
LAM(I)=1000 
PBP(I)=0 
RLC(I)=1 

ELSEIFC  C0G2( I ) .  EQ.  ' 4D ' . OR. C0G2( I ) .  EQ. ' 3D ' . OR.  C0G2 ( I ) .  EQ.  ' 2D ' .  OR. 
*C0G2(I).EQ.  '1D'.0R.C0G2(I).EQ.  '0D')THEN 
RMAX(I)=.  30 
LAM(I)=2000 
PBP(I)=0 
RLC(I)=1 

ELSEIFC C0G2( I ).EQ.  '4E'.0R.  C0G2(I).EQ.  ’3E'.0R.  C0G2(I).EQ.  '2E'.0R. 
*C0G2( I ) .  EQ. ' IE ' .  OR.  C0G2( I ) .  EQ.  ' OE ' )THEN 
RMAX(I)=.  40 
LAM(I)=2000 
PBP(I)=0 


ELSEIF(C0G2(I).EQ.  '4C' .  OR.  C0G2(I).  EQ.  'SC'.OR.  COG2(I).EQ.  '2C')THEN 
RMAX(I)=. 50 
LAM(I)=500 
PBP(I)=0 
RLC(I)=1 

ELSEIFC  C0G2( I ) .  EQ.  ' 1C ' .  OR.  C0G2( I ) .  EQ.  ’ OC ' )THEN 
RMAX(I)=.  99 
LAM(I)=500 
PBP(I)=20 
RLC(I)=0 

ELSEIFC C0G2( I ).EQ.  'AF'.OR.  C0G2(I).EQ.  '3F'.0R.  C0G2(I).EQ.  '2F')THEN 
RMAX(I)=.  50 
LAM(I)=750 
PBP(I)=0 
RLC(I)=1 

ELSEIFC C0G2( I ).EQ. ’ IF' .  OR. C0G2C I).  EQ.  'OF')THEN 
RMAXCI)=.  99 
LAMCI)=500 
PBPCI)=20 
RLCCI)=0 

ELSEIFCC0G2CI).EQ.  'Nl ' .  OR.  C0G2CI).  EQ.  ’N2'.0R.  C0G2CI).EQ.  'N3')THEN 
RMAXCI)=.  35 
LAMCI)=1500 
PBPCI)=0 
RLCCI)=1 
ELSE 

RMAXCI)=.  35 
LAMCI)=2500 
PBPCI)=0 
RLCCI)=1 
ENDIF 
ST0PCI)=0 

IFCRFCD.EQ.  0.  0)RFCI)=DCI) 

IFC RFC  I) .  GT.  DC  I ) )RFC  I )=DC  I ) 

IFCPCLTCD.GT.  18.  )PCLTCI)=18. 

IFC PCLTC I) .  EQ.  0.  0)PCLTC I)=5.  0 
IFCCllCD.GT.  999999.  )C11C I)=999999. 

IFCCllCI).  LT.  .  01)C11CI)=C1CI) 

IFCASICI).  LT. 0.  0)AS1CI)=0.  0 

IFC ASIC  I).  GT.  250000.  )AS1CI)=250000. 

IFCOBSOCD.LE.  .01)0BS0CI)=.  01 
IFCOBSOCI).  GT.  1.  0)0BS0CI)=1.  0 
IFCECI).LT.  .001)ECI)=.  001 
UPLIM=  AMAXOC 100000 , C  NINTC 16*DC I ) ) ) ) 

PPRIYRC I )=PPR1C I)+PPR2C I )+PPR3C I )+PPR4C I ) 

IFCOHRFICI).LT.  0)OHRFICI)=0 

IFCOHRFICI).  GT.  UPLIM)OHRFICI)=UPLIM 

IFC  OHNRFI C I ) .  LT.  0 ) OHNRFI C I )=0 

IFC  OHNRFI C  T ) . GT.  UPLIM) OHNRFI C I )=UPLIM 

IFCADCI).LT. 0)ADCI)=0 

IFCADCI).GT.  UPLIM)ADCI)=UPLIM 

IFCCDCI).LT. 0)CDCI)=0 

IFCCDCI).GT.  UPLIM)CDCI)=UPLIM 

IFCB0CI).LT.0)B0CI)=0 

IFC  BOC I ) .  GT.  UPLIM) BOC I )=UPLIM 

IFCRESVCI).  LT.  0)RESVCI)=0 
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IF(RESV(I).GT.  UPLIM)RZSV(I)=UPLIM 

IF( PPRLT( I ) .  LT.  0 ) PPRLT( I )=0 

IF( PPRLTC I) .  GT.  UPLIM) PPRLT( I )=UPLIM 

IF( PPRIYRC I) .  LT.  0)PPR1YR( I )=0 

IF( PPRIYRC I ) .  GT.  UPLIM) PPRIYRC I )=UPLIM 

ZCI)=DCI)*PCLTCI) 

IFCZCI).  LE.  0.  0)THEN 
ZCI)=0.  0 
PVARCI)=0.  0 
GO  TO  5 
ENDIF 

4  VTMR=150. 

VTMRP=PVARCI)/ZCI) 

IFCVTMRP.  GE.  VTMR)PVARCI)=0.  0 
IFCPVARC • ).NE.  0.  0)G0  TO  5 
PVARCI)=4.  85*ZCI)**1.  5 
5  CONTINUE 
MSRTGG=0.  0 

CALL  ICPM0DCNI,H,B,X,Z,C1,C11,D,QQ,NRPR,RF,PVAR,L0T, 

*DMAD , E , LAM , MARK , PBP , SLC , RMIN , RMAX , RLC , NSO , AS 1 , 

*QQR , QMIN , OBSO , PRIND , HC , POC ) 

VSL=0.  0 
VAT=0. 0 
DO  10  1=1, NI 

*  PRINT*,  QQCD.QQRCI) 

SLCI)=XCI)-QQRCI)-ZCI) 

XUICPCI)=XCI) 

VSL=VSL+C11CI)*SLCI) 

VAT=VAT+C 1 1 C I ) *D  C I ) *4 . 

10  CONTINUE 

SLD=VSL*365.  /VAT 

CALL  PRT0UTC1,NAME1,QM,B,QQ,QRI,QQR,N,NI,NN,X,Z,C11,D,MSRTGG, 
*COGG 1 , C0GG2 , NPO , TOV , QMR , PBP , PVAR , SLD , PCLT , NNN , NN 1 ) 

M0D=1 

CALL  ICPSMACNI, C0GG1,D,RF,DMAD,X,QQ, QMIN, Z, PBP, PVAR, H, RMAX, LAM, 
*RMIN, MARK, SLC, PCLT, Cl, E, MOD, QRI, NSO, 

*NRPR) 

MSRTGG=T0VC2) 

BUDGG=B 
DO  41  1=1, NI 

IFCMQ.EQ.5)QCI)=QQCI) 

IFCMQ.EQ.  5)G0  TO  41 
41  CONTINUE 
M0D=3 

CALL  MODOPTCNI, NNN, B,M0DMST,X,Z,D,QQR,C11, STOP, MSRTGG, MOD, PBP, 
*PVAR,Q,NN1) 

*****  COMPUTE  ROPCD 

DO  7777  1=  1,  NI 
ROPCI)=XCI)-QCI) 

XMSRTCI)=XCI) 

7777  CONTINUE 
VSL=0. 0 
VAT=0. 0 
DO  43  1=1, NI 
SLCI)=XCI)-QCI)-ZCI) 

IFC  SLC I ) .  LT.  0.  0 ) SLC I )=0. 0 
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VSL=VSL+C 1 1 ( I ) *SL( I) 

43  VAT=VAT+C11(I)*D(I)*4. 

SLD=VSL*365.  /VAT 

CALL  PRT0UT(3,NAME3,QM,B,Q,QRI,QQR,N,NI,NN,X,Z,C11,D,MSRTGG, 
*COGG 1 , C0GG2 , NPO , TOV , QMR , PBP , PVAR , SLD , PCLT , NNN , NNl ) 

CALL  ICPSMACNI, COGGl,D,RF,DMAD,X,QQ,QMIN,Z, PBP, PVAR, H.RMAX, LAM, 
*RMIN , MARK , SLC , PCLT , C 1 , E , MOD , QRI , NSO , 

*NRPR) 

CHRNT0T=0.  0 
PRINT* 

WRITE  (6,8184) 

8184  FORMATC O', 'ITEM  CHURN  DMD  UNIT  PRICE  CHURN  PRICE') 

DO  8181  1=1,  NI 

XCHURN(I)  =  XUICP(I)“XMSRT(I) 

CHRNCSTC I )=XCHURN( I )*C 1 1( I ) 

CHRNTOT=CHRNTOT+CHRNCST( I ) 

WRITE(11,8186)XCHURN(I) 

8186  FORMAT  (IX, 18) 

WRITE( 6,8183)1, XCHURN( I),D(I),C11(I), CHRNCST( I ) 

8183  FORMAT( 'o' ,14,18, FIO.  2,F13.  2,F16.  2) 

8181  CONTINUE 

WRITE  (6,8185)CHRNTOT 

8185  FORMAT( ' 0 ' , ' SUM  OF  CHURN  PRICES  IS  $  ' ,F18. 2) 

*1500  CONTINUE 

99  STOP 
END 

*THE  INPUT  FOR  THE  READ  STATEMENT  OF  THE  MAIN  PROGRAM  IS  IN  THE  CARES 
*( COMPUTATION  AND  RESEARCH  EVALUATION  SYSTEM)  III  FORMAT. 

* 


*VARIABLE  NAME 

DATA 

DEN(S) 

CARES  CC# 

*C0G1 

COG 

C003 

1-2 

*C0G2 

COG 

C003W 

3-4 

*MARK 

MARK  CODE 

5 

*SN 

STOCK  NUMBER 

D046D 

6-14 

*SLC 

SHELF  LIFE  CODE 

C028 

16 

*PRIND 

PROCUREMENT  INDICATOR 

D025E 

18 

*0BS0 

OBSOLESCENCE  RATE 

B057 

25  -  27 

*E 

ESSENTIALITY  CODE 

C008C 

28  -  30 

*PCLT 

PROCUREMENT  LEADTIME 

BOllA 

31  -  34 

*NRPR 

NUMBER  OF  POLICY  RECEIVERS 

42  -  45 

*C1 

REPLACEMENT  PRICE 

B055 

46  -  55 

*C11 

UNIT  PRICE 

B053 

56  -  65 

*D 

QUARTERLY  DEMAND 

B074 

66  -  75 

*PVAR 

PROCUREMENT  VARIANCE 

B019A 

86  -  95 

*RF 

REQUISITION  FREQUENCY 

A023B 

96  -  105 

*RESV 

RESERVATIONS 

A013A 

106  -  113 

*B0 

BACKORDERS 

AOll 

114  -  121 

*PPRLT 

PPRS  DURING  LEADTIME 

122  -  126 

*PPR1 

PPRS  1ST  QTR  AFTER  LT 

127  -  131 

*PPR2 

..  2ND  "  "  " 

132  -  136 

*PPR3 

II  gjyj  II  It  II 

137  -  141 

*PPR4 

"  4TH  "  "  " 

142  -  146 

*0HRFI 

ON -HAND  SYSTEM  RFI 

147  -  154 

*0HNRFI 

ON-HAND  SYSTEM  NRFI 

155  -  162 

*AD 

AWARDED  DUES 

163  -  170 
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*CD 

COMMITTED  DUES 

171  -  178 

*AS1 

SET-UP  COST 

BOSS 

179  -  186 

*DMAD 

DEMAND  MAD  A019(2)+A019A(2) 

225  -  234 

*LOT 

LIFE  OF  TYPE  QTY 

B070 

235  -  242 

*NSO 

REORDER  LVL  LOW  LIMIT  QTY 

B020 

243  -  250 

*QMIN 

MINIMUM  PRODUCTION  QTY 

B061B 

251  -  255 
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APPENDIX  B.  SUBROUTINE  ICPMOD 


C  ***  THE  CURRENT  I CP  CONSUMABLES  MODEL 

c  ***  this  subroutine  computes  only  the  procurement  q  and  ROP. 

SUBROUTINE  ICPMOD(N,H,B ,X,Z,C1 ,C11 ,D,Q,NRPR,RF,PVAR,LOT, 

*DMAD , E , LAM , MARK , PBP , SLC , RMIN , RMAX , RLC , NSO , AS 1 , QQR , QMIN , 
*OBSO,PRIND,HC,POC) 

INTEGER  N,NRPR(N),LOT(N),MARK(N),Y(1500),R121,Q8,Q9,QMIN(N) 
INTEGER  X(N) ,NSO(N) ,RLC(N) ,PBP(N) ,ROP( 1500) ,R0P1, I, J,K 
CHARACTER*!  PRIND(N) 

REAL  Z(N) ,C1(N) ,D(N) ,RF(N) ,PVAR(N) ,DMAD(N) ,AS1(N) 

REAL  SL( 1500) ,RMIN(N) ,RMAX(N) ,C3( 1500) ,RISK( 1500) ,E(N) ,SLC(N) ,H 
REAL  Q1,Q(N),LAM(N),C11(N) 

REAL  QQR(N) ,AA,AC,RRCT,HC(N) ,OBSO(N) ,POC(N) 

REAL  TEMPI , TEMP2 , T1 , QIA , QIB , QIC , QID , QIE , QISQRT , QIMAX 
REAL  OLP1P2,TP1P2,Z1B23H,B19AOC,TSTRSK,RISKTB(50),TVALU,RISKJ 
REAL  RESLT , VTMTR , PVTMR , QRATIO , VRATIO , VLOG , PROB , SMPROB , B 19PT3 
REAL  B19PT1,B19PT2,B19PT4,B19MIN,B19MAX,B21MIN,B21PT1,B21PT2 
REAL*8  B 
B=0.  0 

*****  RISKTABLE  for  THE  NORMAL  DISTRIBUTION 
RISKTB(l)  =  0.46017:20 
RISKTB(2)  =  0.42074060 
RISKTB(3)  =  0.38206860 
RISKTB(4)  =  0.33457830 
RISKTB(5)  =0.30853750 
RISKTB(6)  =0.27485310 
RISKTB(7)  =  0.24196370 
RISKTB(8)  =  0.21185540 
RISKTB(9)  =0.18406010 
RISKTB(IO)  =  0.15865530 
RISKTB(ll)  =  0. 13566610 
RISKTB(12)  =  0. 11506970 
RISKTB(13)  =  0.09680050 
RISKTB(14)  =  0.08075670 
RISKTB(15)  =  0.06680720 
RISKTB(16)  =  0.05479930 
RISKTB(17)  =  0.04456550 
RISKTB(18)  =  0.03593030 
RISKTB(19)  =  0.02871660 
RISKTB(20)  =  0.02275010 
RISKTB(21)  =  0.01786440 
RISKTB(22)  =  0.01390340 
RISKTB(23)  =  0.01072410 
RISKTB(24)  =  0.00819750 
RISKTB(25)  =  0.00620970 
RISKTB(26)  =  0.00466120 
RISKTB(27)  =  0.00346700 
RISKTB(28)  =  0.00255510 
RISKTB(29)  =  0.00186580 
RISKTB(30)  =  0.00134990 
RISKTB(31)  =  0.00096760 
RISKTB(32)  =  0.00068710 
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RISKTB(33)  =  0.00048340 
RISKTB(34)  =  0.00033690 
RISKTB(35)  =  0.00023260 
RISKTB(36)  =  0.00015910 
RISKTB(37)  =  0.00010780 
RISKTB(38)  =  0.00007230 
RISKTB(39)  =  0.00004810 
RISKTB(40)  =  0.00003170 
RISKTB(41)  =  0.00002070 
RISKTB(42)  =  0.00001330 
RISKTB(43)  =  0.00000850 
RISKTB(44)  =  0.00000540 
RISKTB(45)  =  0.00000340 
RISKTB(46)  =  0.00000211 
RISKTB(47)  =  0.00000130 
RISKTB(48)  =  0.00000079 
RISKTB(49)  =  0.00000048 
RISKTB(50)  =  0.00000029 
DO  6  1=1, N 

******  COMPUTE  HOLDING  COST  (HC).  HC  =  STORAGE  RATE(.Ol)  +  TIME 
******  preference  RATE(.IO)  +  OBSOLESCENSE  RATE  (OBSO(I)) 
HC(I)=OBSO(I)+.  11 

******  COMPUTE  AND  CONSTRAIN  THE  RISK 
TEMP1=HC(I)*D(I) 

TEMP2=RF( I )*LAM( I)*E( I ) 

IF(TEMP2.EQ.  0.0) THEN 
T1=0. 0 
ELSE 

T1=AMIN1(TEMP1*C1(I)/TEMP2, 999999.0) 

ENDIF 

RISK(I)=T1/(T1+1.0) 

RISK(I)=AMAX1(RISK(I),RMIN(I)) 

RISK( I )=AMIN1(RISK( I ) ,RMAX( I) ) 

******  COMPUTE  PROCUREMENT  ORDER  COST  (POC) 

IF( MARK( I ) .  EQ.  0.  OR.  MARK( I ) .  EQ.  1.  OR.  MARK( I ) .  EQ.  2 )THEN 
P0C(I)=660. 

ELSEIF(MARK( I).  EQ.  3.  OR.  MARK(I).  EQ.  4. AND. (8000**2*HC( I)* 
*(AMAXl((l-RISK(I)),(l-.5)))/8*(660.+ASl(I))).GT.  (C1(I)*D(I)))THEN 
POC(I)=660. 

ELSEIF(PRIND(I).EQ.  'O' . OR.  PRIND( I).  EQ.  'B'.  OR.  PRIND(I).  EQ.  '2') 

*  THEN 

P0C(I)=1940. 

ELSE 

P0C(I)=1970. 

ENDIF 

COMPUTE  AND  CONSTRAIN  REORDER  POINT  (ROP) 

*****  BASIC  REORDER  LEVEL  ***** 

**  THE  TEMPORARY  VARIABLE  NAMES  USED  IN  THE  REMAINDER  OF  THIS 
**  SUBPROGRAM  WERE  USED  IN  ORDER  TO  MATCH  AS  CLOSELY  AS  POSSIBLE 
**  THE  VARIABLE  NAMES  USED  IN  THE  ACTUAL  UICP  PROGRAM 
OLP1P2=1.0-RISK(I) 

TP1P2=RISK(I) 

Z1B23H=Z(I) 

B19AOC=PVAR(I) 

IF(MARK(I).EQ.  0)GO  TO  1500 
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IF  (Z1B23H. LT. PBP(I))GO  TO  1000 
*****  NORMAL  DISTRIBUTION 
500  TSTRSK=TP1P2 

IF(TP1P2.  LE.  0.  5)GO  TO  600 
TSTRSK=1.  0-TP1P2 
600  DO  650  K  =  1,  50 

IF(TSTRSK.  LT.RISKTB(K))GO  TO  650 
TVALU=1.  0*K 

IF  (TSTRSK.EQ.RISKTB(K))GO  TO  700 

J=K-1 

RISKJ=0. 5 

IF( J.  NE.  0)RISKJ=RISKTB( J) 

TVALU=K-( (TSTRSK-RISKTB(K) )/(RISKJ-RISKTB(K) ) ) 
GO  TO  700 
650  CONTINUE 

TVALU=50.  0 

700  IF(TP1P2.  GT.  0.  5)TVALU=-TVALU 

TVALU=0.  1*TVALU 
800  K=0 

RESLT=Z1B23H 

IF(B19A0C.  LE.  0.  0)GO  TO  2500 

IF(TP1P2. NE. 0. 5)RESLT=Z1B23H+(TVALU*SQRT(B19A0C)) 
GO  TO  2500 

*****  negative  BINOMIAL  DISTRIBUTION  ***** 

1000  VTMR=1. 01 

IF(Z1B23H.  NE.  0.  0)VTMR=AMAX3(VTMR,(B19AOC/Z1B23H)) 

PVTMR=VTMR-1.  0 

QRATIO=PVTMR/VTMR 

VRATI0=Z1B23H/PVTMR 

VLOG=VRATIO*ALOG( VTMR) 

IFCVLOG.GT.  6.  9)G0  TO  500 
1100  K=0 

PROB=EXP( -VLOG) 

SMPROB=PROB 
GO  TO  1300 
1200  K=K+1 

PROB=( ( VRATIO+( K- 1 ) ) /K)*QRATIO*PROB 
SMPROB=SMPROB+PROB 
1300  IFCPROB.  LE.  0.  0001 )GO  TO  500 

IF(SMPROB. LT.0LP1P2)G0  TO  1200 
RESLT=1.  0*K 
GO  TO  2500 

*****  POISSON  DISTRIBUTION  ***** 

1500  IF  (Z1B23H.  GT.  6.  9)G0  TO  500 

1510  K=0 

PR0B=EXP(-Z1B23H) 

SMPROB=PROB 
GO  TO  1540 
1520  K=K+1 

PR0B=(Z1B23H/K)*PR0B 
SMPROB=SMPROB+PROB 
1540  IF(PROB.LE.  0.  OOODGO  TO  500 

IF(SMPROB.  LT.0LP1P2)G0  TO  1520 
1550  RESLT=1. 0*K 

GO  TO  2500 
2500  RESLT=RESLT 
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IF  (D(I).GT.  0.  0)GO  TO  3860 
RESLT=AMAX1( Z( I ) ,FLOAT(NSO( I) ) ,0.0) 

CALL  SHFLIF(SL(I),SLC(I)) 

GO  TO  3990 

3860  CALL  SHFLIF(SL(I) ,SLC(I)) 

B19PT3=SL(  I  )-"D(  I  )+Z(  I ) 

IF  (NSO(I).GT. 100000.  0)NSO(I)=0.0 
B19PT1=Z(I)+33.*D(I) 
B19PT2=AMAX1(RESLT,FL0AT(NRPR( I) ) ) 
IF((SL(I)/4.  ).LT.  7.0)GO  TO  3900 
B19PT3=99999999.0 
GO  TO  3920 

3900  B19PT3=B19PT3*D(I) 

3920  B19MIN=AMIN1(B19PT1,B19PT2,B19PT3) 

B19PT4=1.  0*Z(I) 

B19MAX=AMAX1(B19MIN,B19PT4,FL0AT(NS0(I)),0.  0) 
3940  RESLT=B19MAX+0. 999 

3990  R0P(I)=AINT(RESLT) 

6  CONTINUE 

*****  COMPUTE  Q  NEXT 
DO  12  1=1, N 

Q1A=12.  0*D(I) 

Q1B=D(I) 

Q1C=8.  0*(POC(I)+AS1(I))*D(I) 

Q1D=HC( I )*C 1( I )*AMAX1( ( 1 -RISK( I ) ) , .  5 ) 

Q1E=0.  0 

IFCQID.GT.  0.  0)THEN 
Q1E=Q1C/Q1D 
ENDIF 

IF(Q1E.LE.0.0)THEN 
Q1SQRT=0. 0 
ELSE 

Q1SQRT=SQRT(Q1E) 

ENDIF 

Q1MAX=AMAX1(Q1B ,Q1SQRT) 

Q1=AMIN1(Q1MAX,Q1A) 

Ql=Ql+0.  999 

IF((SL(I)/4.  ).LT.  7.0)G0  TO  4020 
Q(I)=Ql+0.  999 
GO  TO  4100 

4020  B21PT1=AMAX1(R0P(I)-Z(I),0. 0) 

B21PT2=SL( I)*D( I) -B21PT1 
B21MIN=AMIN1(B21PT2,Q1) 

Q(I)=B21MIN+0.  999 
*4100  IF(Q(I).LT.  1.0)Q(I)=1.0 

4100  Q(I)=AINT(Q(I)) 

QQR(I)=Q(I) 

X( I )=ROP( I )+IFIX( QQR( I )+0.  999) 

B=B+C 11(1) *FLOAT( X( I ) ) 

12  CONTINUE 
RETURN 
END 

* 

*VARIABLES  NOT  IDENTIFIED  IN  APPENDIX  1 
* 

*VARIABLE  NAME  DATA 
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*N 

*B 

*X 

*Q 

*LAM 

*PBP 

*RMIN 

*RMAX 

*RLC 

*QQR 

*HC 


NUMBER  OF  ITEMS  BEING  PROCESSED 
(NOT  USED) 

BUDGET  (INVESTMENT  LEVEL) 

DEPTH 

ECONOMIC  ORDER  QUANTITY 

LAMBDA  (GIVEN  IN  THE  MAIN  PROGRAM) 

PROBABILITY  BREAK  POINT  " 

MINIMUM  RISK  " 

MAXIMUM  RISK  " 

REORDER  LEVEL  CONSTRAINT  ” 

ECONOMIC  ORDER  QUANTITY 
HOLDING  COST 
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APPENDIX  C.  SUBROUTINES/FUNCTIONS  -  MODOPT,  MODMST 

AND  GOALM 

*  IN  THIS  APPENDIX  THREE  SUBPROGRAMS  HAVE  BEEN  LISTED  TOGETHER 
♦BECAUSE  THEY  FUNCTION  AS  ONE  SUBPROGRAM  AND  COULD  HAVE  BEEN  WRITTEN 
♦AS  SUCH. 


SUBROUTINE  MODOPT( N , NNN , B , AMODEL , X , ZN , D , QR , C 1 1 , STOP , GOALG , MOD , PBP , 
♦PVAR.Q.NNl) 

C  — -  PERFORM  OPTIMAL  ALLOCATION  FOR  GIVEN  MODEL  USING 
C  ---  MARGINAL  ANALYSIS  METHOD  AND  LOWER  BOUNDING. 

C  ---  N=NO.  ITEMS 

C  ---  B=INVESTMENT  LEVEL  OF  STOCK  FUND 

C  ---  AMODEL=ENTRY  POINT  FOR  MODEL  TO  USE  (STANDARDIZED  ARGUMENTS) 

C  ---  X=OPTIMAL  ALLOCATIONS  PER  ITEM 
C  ---  ZN=  MEAN  DEMAND  DURING  RESUPPLY  TIME  OR  PPV 
C  ...  Cl=PROCUREMENT  COST  FOR  EACH  ITEM 
C  ---  RR=WORK  VECTOR  TO  STORE  RATIOS 

INTEGER  N, NNN, K,MK, STEP, X(N) ,STOP(N) ,PBP(N) , INDEX( 1500) ,XL(1500) 
REAL  ZN(N) ,QR(N) ,D(N) ,MR,RR( 1500) ,SR,TRY,PVAR(N) ,Q(N) 

REAL  MODNSF, MODMST, Cll(N), GOALG 
REAL*8  B.BR,BUDLIM 
1=1 
SR=0. 

KK=1 

NN1=0 

C  ---  INITIALIZE 
BUDLIM=B 
B=0. 

DO  2  1=1, N 

C  ---  INITIALIZE  STOP  BEFORE  OPTIMIZING  ON  INVESTMENT  LEVEL(STOP=l  MEANS 
C  ---  THAT  THE  LEVEL  HAS  HIT  THE  ITEM  MSRT  BOUND). 

******XHESE  FIVE  STEPS  PUT  X  AT  ROP**^**^ 

♦  RP(I)  =  X(I)-IFIX(QR(I)) 

♦  IF  (X(I).GT.  1)  THEN 

♦  X(I)  =  RP(I) 

♦  ENDIF 

♦  B  =  B+(X(I)^C11(I)) 

♦^^♦♦♦♦THESE  FOUR  STEPS  PUT  X  AT  1  OR 

♦  IF(X(I).NE.  0)THEN 

♦  X(I)=1 

♦  B=B+C11(I) 

♦  ENDIF 

♦♦♦♦♦♦this  step  puts  x  at  zero 
x(I)=0 
ST0P(I)=0 

C  INDEX(I)=0 

RR(I)=AMODEL(ZN(I),D(I),QR(I),Cll(I),X(I)+l,STOP(I),PBP(I), 

♦  PVAR(I)) 

2  CONTINUE 
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no  n  nooo 


12  STEP=0 

C  —  DO  UNTIL  MSRT  OR  BUDGET  GOAL  IS  REACHED 
20  CONTINUE 
STEP=STEP+1 
MK=0 

MR=1S000000. 

IF(M0D. EQ. 3)MR=0. 

DO  30  K=1,N 

IF(STOP(K).EQ.  1)G0  TO  30 
*  IF(INDEX(I).EQ.  1)G0  TO  30 

IF(RR(K).GE.MR)GO  TO  30 
MR=RR(K) 

MK=K 
30  CONTINUE 

IF(MK  .  EQ.  0)  GO  TO  50 

C  ---  ALLOCATE  ONE  MORE  UNIT  OF  ITEM  MK  IF  POSSIBLE. 

B=B+C11(MK) 

IF(B. GE. BUDLIM)THEN 
B=B-C11(MK) 

GO  TO  50 
ENDIF 

X(MK)=X(MK)+1 

---  NEXT  CHECK  TO  SEE  IF  GOAL  HAS  BEEN  ATTAINED. 

CALL  GOALM( X , N . ZN , D , QR , TRY , GOALG , PBP , PVAR ) 

IFCTRY.  EQ.  0.  0)GO  TO  50 
SR=MR 

RR(MK)=AMODEL(ZN(MK) ,D(MK) ,QR(MK) ,C11(MK) ,X(MK)+l,STOP(MK) , 
*PBP(MK),PVAR(MK)) 

GO  TO  20 
50  RETURN 
END 

---  ROUTINE  TO  MINIMIZE  MEAN  SUPPLY  RESPONSE  TIME  —  MSRT 
REAL  FUNCTION  MODMST(ZZ,D,QR,C,K, STOP, PBP, PVAR) 

COMPUTE  MARGIN  ANALYSIS  RATIOS  ASSUMING 
POISSON  DEMAND. 

REAL  ZZ,QR,C, MSRT, MSRTD,D, PVAR 
INTEGER  K, STOP, PBP 
IFCZZ.  NE.  0. 0)G0  TO  10 
M0DMST=0.  0 
MSRT=-0.  0001 
GO  TO  12 

10  TS1=TWUS(ZZ,QR,K-1,PBP,PVAR,MARK) 

TS2=TWUS( ZZ , QR , K , PBP , PVAR , MARK) 

M0DMST=(TS2-TS1)/C 
—  NOTE.  MODMST  WILL  BE  NEGATIVE 

MSRT=TVUS( ZZ , QR , K , PBP , PVAR, MARK)/D 
12  MSRTD=91.  *MSRT 

IF(MSRTD.  LT.  0.  01)STOP=1 
RETURN 
END 

— -  ROUTINE  TO  SEE  IF  GOAL  HAS  BEEN  ATTAINED. 

SUBROUTINE  GOALM( X , N , Z , D , QR , TRY , MSRTG , PBP , PVAR ) 

INTEGER  N,X(N),XI,PBP(N) 

REAL  Z(N) ,VO( 1500) ,TV,SLT,D(N) ,QR(N) ,MSRTG,TRY 
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REAL  MSRTC,MSRT,PVAR(N) 

TV=0. 0 
SLT=0. 

DO  10  1=1, N 

IF(Z(I).EQ.  0.  0)GO  TO  16 
XI=X(I) 

SLT  =  SLT  +  D(I) 

14  MSRTC=TWUS( Z ( I ) , QR( I ) ,XI , PBP( I ) , PVAR( I ) .MARK) /D( I ) 

MSRT=AMAX1(MSRTC,0.  0) 

GO  TO  18 
16  MSRT=0. 0 

18  VO(I)  =  91.*MSRT 

TV  =  TV  +  VO(I)*D(I) 

10  CONTINUE 

TV=TV/SLT 
TRY  =  1.  0 

IF(TV.  LE.MSRTG)TRY  =  0.0 
22  RETURN 
END 

* 

*VARIABLES  NOT  DEFINED  IN  PREVIOUS  APPENDICES 
* 

*VARIABLE  NAME  DATA 

(NOT  USED) 

MODMST 

LEADTIME  DEMAND 
ECONOMIC  ORDER  QUANTITY 

1  OR  0  (INDICATES  THAT  THE  ITEM  HAS  MET  THE  MSRT  GOAL) 
MSRT  GOAL  (EQUALS  THE  MSRT  THE  UICP  MODEL  COMPUTED  FOR 
THAT  FOUR  DIGIT  COG 
(NOT  USED) 

(NOT  USED) 

LEADTIME  DEMAND 
UNIT  PRICE 
INTERMEDIATE  DEPTH 

1  OR  0  (INDICATES  THAT  THE  FOUR  DIGIT  COG  HAS  MET  THE 
MSRT  GOAL) 


*NNN 

*AMODEL 

*ZN 

*QR 

*STOP 

*GOALG 

* 

*MOD 

*NN1 

*ZZ 

*C 

*K 

*TRY 

* 
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APPENDIX  D.  SUBROUTINE  SHFLIF 


SUBROUTINE  SHFLIF(SL,SLC) 

REAL*4  A/'A  7.C/’C  ’/,D/'D  '/,E/’E  '/, 

*F/'F  7, 

*G/'G  7.h/’h  7,J/’J  7,k/’k  7,l/'l  7.m/'m  7. 

*N/'N  7,P/’P  7,Q/'Q  7,R/’R  7,x/'x  7. s/'s  '/, 

*Al/'l  '/,A2/’2  '/,A3/’3  */,A4/'4  '/,A5/'5  ’/, 

*A6/'6  7,A7/7  7.A8/’8  '/,A9/'9  '/,A0/'0  '/ 

REAL*4  SL.SLC 
IF(SLC.EQ.  A0)SL=100. 

IF(SLC.EQ.  AO)GO  TO  50 
IF(SLC.  EQ.A)SL=l./3. 

IF(SLC.EQ.  A)GO  TO  50 
IF(SLC.  EQ.  B)SL=2. /3. 

IF(SLC.  EQ.  B)GO  TO  50 

IF(SLC.  EQ.  C.  OR.  SLC.  EQ.  A1)SL=3.  /3. 

IF(SLC.EQ.  C.  OR.  SLC.EQ.  Al)GO  TO  50 
IF(SLC.EQ.D)SL=4.  /3. 

IF(SLC.EQ.  D)GO  TO  50 
IF(SLC.  EQ.  E)SI,=5.  /3. 

IF(SLC.  EQ.  E)GO  TO  50 

IF(SLC.  EQ.  F.  OR.  SLC.  EQ.  A2)SL=6.  /3. 

IF(SLC.  EQ.F.  OR.  SLC.EQ.  A2)G0  TO  50 
IF(SLC.  EQ.  G.  OR.  SLC.  EQ.  A3)SL=9.  /3. 

IF(SLC.  EQ.G.  OR.  SLC.EQ.  A3)G0  TO  50 
IF(SLC.  EQ.  H.  OR.  SLC.  EQ.  A4)SL=12.  /3. 

IF(  SLC.  EQ.  H.  OR.  SLC.  EQ.  A4)G0  TO  50 
IF(SLC.EQ.  J)SL=15.  /3. 

IF(SLC.  EQ.  J)GO  TO  50 

IF( SLC.  EQ.  K.  OR.  SLC.  EQ.  A5)SL=18.  /3. 

IFCSLC.EQ.K.  OR.  SLC.EQ.  A5)G0  TO  50 
IF(SLC.EQ.  L)SL=21.  /3. 

IF(SLC.  EQ.  L)GO  TO  50 

IF(  SLC.  EQ.  M.  OR.  SLC.  EQ.  A6)SL==24.  /3. 

IF( SLC.  EQ.  M.  OR.  SLC.  EQ. A6)G0  TO  50 
IF(SLC.  EQ.  N)SL=27.  /3. 

IFCSLC.EQ.  N)GO  TO  50 
IF(SLC.EQ.  P)SL=30.  /3. 

IFCSLC.EQ.  P)GO  TO  50 

IFCSLC.  EQ.  Q.  OR.  SLC.  EQ.  A7)SL=36.  /3. 

IFCSLC.EQ.Q.  OR.  SLC.EQ.  A7)G0  TO  50 
IF( SLC.  EQ,  R.  OR.  SLC.  EQ. A8)SLs48.  /3. 

IFCSLC.  EQ.R,  OR.  SLC.EQ.  A8)G0  TO  50 
IFCSLC.  EQ.  S.  OR.  SLC,  EQ.  A9)SL=60.  /3. 

IFCSLC.EQ.  S.  OR.  SLC.  EQ.A9)G0  TO  50 
SL=100. 

50  RETURN 
END 

* 

★VARIABLE  NOT  DEFINED  IN  PREVIOUS  APPENDICES 
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*VARIABLE  NAME  DATA 

*SL  ITEM  SHELF  LIFE 


APPENDIX  E.  REAL  FUNCTION  TWUS 

C  ---  ROUTINE  TO  CALCULATE  THE  EXP  TIME  VTID  UNITS  SHORT  FOR  K  UNITS 
REAL  FUNCTION  TWUS(ZZ,QR,K,PBP,PVAR,MARK) 

REAL  ZZ,P1,P2,SW,RP,QR,P3,P4.P5,P6 

REAL  CCD 1 , CCD2 , CCD3 , CCD4 , CCDS , CCD6 , BETAl , BETA2 

REAL*8  PHI 1 , PHI 2 , DCD 1 , DCD2 , D 1 , D2 , T1 , T2 , Z , DBETAl , DBETA2 , S , R , PV 

REAL*8  CD1,CD2,CD3,CD4,CD5.CD6 

INTEGER  K,MARK,PBP 

SW=FLOAT(K) 

S=SW 

PV=PVAR 

Z=ZZ 

RP=SW-QR 

R=RP 

KRP=K-IFIX(QR+.  5) 

IF(ZZ.GE.FLOAT(PBP))GO  TO  20 
CALL  CDFP(ZZ,KRP-1,P1,CD1,CCD1) 

CALL  CDFP(ZZ,KRP,P2,CD2,CCD2) 

CALL  CDFP  (ZZ,KRP+1,P3,CD3,CCD3) 

CALL  CDFP(ZZ,K-1,P4,CD4,CCD4) 

CALL  CDFP(ZZ,K,P5,CD5,CCD5) 

CALL  CDFP(ZZ,K+1,P6,CD6,CCD6) 

IF  (CCDl.LT.  0.000001)  GO  TO  10 

BETA 1=( CCD1*ZZ**2 ) /2.  -CCD2*ZZ*KRP+CCD3*KRP*(KRP+1 ) /2. 

IF(BETA1.  LT.  0.  00000 1)BETA1=0.  0 
IF  (CCD4.GE.  0.000001)  GO  TO  7 
BETA2=0.  0 
GO  TO  8 

7  BETA2=(CCD4*ZZ**2)/2.  -CCD5*ZZ*K+CCD6*K*(K+l)/2. 

IF(BETA2. LT. 0. 000001)BETA2=0.  0 

8  TW'US=(  BETAl  "BETA2)/QR 
RETURN 

10  ■IVUS=0.  0 
RETURN 

20  T1=(R-Z)/DSQRT(PV) 

IF(DABS(T1).LE.  7.  )GO  TO  21 
PHI 1=0. 0 

IF(T1.  GT.  7.  )D1=0.  0 
IF(T1.LT.  (-7.0))D1=1.0 
GO  TO  22 

21  PHIl=(DEXP(-(Tl**2)/2.  ) )/SQRT(2.  *3.  14159265) 

CALL  HDN0RDCT1,CD1) 

Dl=l.  0-CDl 

*  PRINT*,  PV,D1,T1 

*  PRINT*,  PHIl 

22  DBETA1=PV*(D1*(1.  0+Tl**2)-Tl*PHIl)/2. 

T2=(S-Z)/DSQRT(PV) 

IF(DABS(T2).LE.  7.  )GO  TO  23 

PHI 2=0.  0 

IF(T2.GT.  7.0)D2=0.  0 
IF(T2.LT.  (-7.0))D2=1.0 
GO  TO  25 
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23  PHI2=(DEXP(-(T2**2)/2.  ))/SQRT(2.  *3.  14159265) 
CALL  MDN0RD(T2,CD2) 

D2=l.  0-CD2 

25  DBETA2=PV*(D2*(1.  0+T2**2) -T2*PHI2)/2. 

TWUS=( DBETAl -DBETA2 ) /QR 

RETURN 

END 

* 

*ALL  VARIABLES  DEFINED  IN  PREVIOUS  APPENDICES 
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APPENDIX  F.  SUBROUTINE  ICPSMA 

SUBROUTINE  ICPSMA(N,C0GG1,D,RF,DMAD,X,QQ,QMIN,Z,PBP,PVAR,H,RMAX, 
*LAM , RMIN , MARK , SLC , PCLT , C 1 , E , 

*MOD,QR,NSO,NRPR) 

INTEGER  DEF , DEFR , ROP , ROPA , QA , QMIN( N ) , XE , ZI , Q2 , X23 , X24 
INTEGER  XH,X2,X21,B019B,MARK(N),QS,Q21,Q2C,Y( 1500), X(N), MOD 
INTEGER  NSO(N) ,NRPR(N) ,PBP(N) ,ROPR( 1500) 

CHARACTER*4  COGGl(N) 

REAL  Z(N) ,DMAD(N) ,D(N) ,PVAR(N) ,PP( 1500) 

REAL  TVEBC 1500) ,ETN( 1500) ,RF(N) ,PCLT(N) 

REAL  H 

REAL  E(N) ,C1(N) ,SLC(N) ,QR(N) ,QQ(N) ,EF,EG 
REAL  Size  1500), RMAX(N),RMIN(N),LAM(N) 

REAL*8  ZZ,PV,R,T1,PHI1,P1,CPHI1,USHRT1,WEB1,Q,T2,PHI2,P2,CPHI2 
REAL*8  USHRT2 ,WEB2 ,TR1 ,TR2 ,PHIR1 ,PHIR2 ,CPHIR1 ,CPHIR2 ,PR1 ,PR2 
REAL*8  EUS,ERUS3,ERUS4,REPQ 
DO  100  1=1, N 

C  ***  COMPUTE  REQUISITIONS  SHORT  FOR  PROCUREMENT  CYCLE,  ETN 
ZZ=Z(I) 

PV=PVAR(I) 

ROP=X( I ) - IFIX( QQ( I )+0.  999 ) 

IF  (D(I).  EQ.  0.  0)THEN 
ETN(I)=0.  0 
TWEB(I)  =0.0 
GO  TO  100 
ENDIF 

3  IF(DMAD(I).NE.  0.0)G0  TO  5 
DMADC I )=( 1.  39*D( I )**0.  75 )**2 
GO  TO  5 

5  EF=1.  57*DMAD(I)/(2.*D(I)) 

W1=0. 154 
EG=W1*D(I)*0.  5 
DEF=IFIX(EF+EG+0.  999) 

ROPA=ROP-DEF 
QA=IFIX(QQ(I)+0.  5)+DEF 
XE=MAXO(QA,QMIN(I)) 

Q=FLOAT(XE) 

IFCROPA. GT. 0)GO  TO  10 
ZI=IFIX(Z(I)+0.  999) 

KN=MINO(ZI,XE) 

ETN( I )=( 4. *RF( I )*FLOAT(KN) )/FLOAT( XE) 

TVEBC I )=ETN( I )*( PCLTC I )+0.  5*W1) /4. 

GO  TO  100 

10  IF(Z(I).LT.FLOAT(PBP(I)))GO  TO  20 
R=FLOAT(ROPA) 

T1=(R-ZZ)/DSQRT(PV) 

IF(DABS(T1).GT.  7.  )GO  TO  11 
PHIl=(DEXP(-(Tl**2)/2. ))/SQRT(2.*3.  14159265) 

CALL  MDN0RD(T1,P1) 

CPHI  1=1.0  -  PI 
GO  TO  12 

11  IFCTl.GT.  7.  )CPHI1=0.0 
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IF(T1.LT.  (-7.  ))CPHI1=1.  0 
PHI 1=0.  0 

12  USHRT1=( DSQRT( PV) )*( PHI 1 -T1*CPHI 1 ) 

WEB 1=PV*( CPHI 1*( 1.  +T1**2 ) -T1*PHI 1 ) 

T2=(R+Q-ZZ)/DSQRTCPV) 

AB=(R+Q) 

IF(DABS(T2).GT.  7.  )G0  TO  13 
PHI2=(DEXP(-(T2**2)/2. ))/SQRT(2.*3.  14159265) 

CALL  MDN0RD(T2,P2) 

CPHI2=1.  0  -  P2 
GO  TO  14 

13  IF(T2.  GT.  7.  )CPHI2=0.  0 
IF(T2.LT.  (-7.  ))CPHI2=1.0 
PHI 2=0.  0 

14  USHRT2=(DSQRT(PV))*(PHI2-T2*CPHI2) 

WEB2=PV*(CPHI2*( 1.  +T2**2)-T2*PHI2) 

GO  TO  50 

20  BQ=PVAR(I)/Z(I) 

IF(BQ.  GE.  1.  0.  AND.  BQ.  LT.  1.  000001)BQ=1.  000001 
BK=Z(I)/(BQ-1.  0) 

PP(1)=1.  O/CBO^'-BK) 

USHRT1=Z( I ) -FLOAT( ROPA)*( 1.  -PP( 1) ) 

WEB1=Z( I )**2+PVAR( I ) -2.  *Z( I)*FLOAT(ROPA)+( (FLOAT(ROPA) )**2)*( 1.  - 
*PP(1)) 

IFCROPA.  LE.  l)GO  TO  23 

KM1=R0PA-1 

DO  22  J=1,KM1 

PP(J+1)=((BK+J-1)/(J))*((BQ-1.  )/BQ)*PP(J) 

USHRT1=USHRT1 -FLOAT( J-ROPA)*PP( J+1 ) 

22  WEB 1=WEB 1 - ( FLOAT( J-R0PA)**2 )*PP( J+1 ) 

23  KM2=R0PA+XE-1 
IF(KM2.LE.  l)GO  TO  50 

USHRT2=Z( I) -FLOAT(KM2+l )*( 1.  -PP( 1) ) 

WEB2=Z( I)**2+PVAR( I) -2.  *Z( I)*FLOAT(KM2+l)+( (FLOAT(KM2+l) )**2)*( 1. 
*PP(1)) 

DO  24  J=1,KM2 

PP(J+1)=((BK+J-1)/(J))*((BQ-1.  )/BQ)*PP(J) 

USHRT2=USHRT2 - ( FLOAT( J -KM2 - 1) )*PP( J+ 1 ) 

24  WEB2=WEB2 -( ( FLOAT( J-KM2- 1 ) )**2)*PP( J+1) 

50  EUS=DMIN1((USHRT1-USHRT2),Q) 

ETN(I)=EUS*4.  *RF(I)/Q 
EW=( WEB 1 -WEB2 )*RF( I ) /D( I ) 

TWEB(I)=EW/(2.*Q) 

100  CONTINUE 

C  ***  COMPUTE  AGGREGATE  SMA 

101  SUMET=0.  0 
SUMTWB=0.  0 
SUMRF=0.  0 

DO  102  1=1, N 
SUMET=SUMET+ETN(I) 

SUMTWB=SUMTWB+TWEB ( I ) 

102  SUMRF=SUMRF+RF(I) 

SMA2=1 - ( SUMET) / ( 4.  *SUMRF) 

C  ***  COMPUTE  AVERAGE  DAYS  DELAY 

ADDA2=365.  *(SUMTWB)/(4.  +SUMRF) 

C  ***  COMPUTE  ADD  FOR  BACKORDERS 
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ADDB02=ADDA2/(1.  -SMA2) 

SMA2=SMA2*100.  0 
PRINT* 

WRITE(6,149) 

149  FORMAT(*0’ ,3X,' CARES  REQUISITION-BASED  PERFORMANCE  MEASURES:') 
WRITE ( 6 , 15 0 ) SMA2 , ADDA2 , ADDB02 

150  FORMATC'O’ .5X,'SMA=' ,F8.  4,5X, ’ ADD=’ ,F8.  4,5X, ' ADDBO=' ,F10.4) 

DO  170  1=1, N 

IF(RF(I).  EQ.  0.  0)RF(I)=1.  0 
SIZ(I)=D(I)/RF(I) 

ETN(I)=ETN(I)*SIZ(I) 

TWEB( I )=TWEB( I )*SIZ( I ) 

170  CONTINUE 
SUMET=0.  0 
SUMTWB=0.  0 
SUMRD=0.  0 
DO  180  1=1, N 
SUMET=SUMET+ETN( I) 

SUMTWB=SUMTWB+TWEB( I ) 

180  SUMRD=SUMRD+D( I ) 

SMA2=l-(SUMET)/(4.  *SUMRD) 

C  ***  COMPUTE  AVERAGE  DAYS  DELAY 

ADDA2=365.  *(SUMTWB)/(4.  *SUMRD) 

C  ***  COMPUTE  ADD  FOR  BACKORDERS 
ADDB02=ADDA2/(1.  -SMA7) 

PRINT* 

SMA2=SMA2*100.  0 
WRITE(6,190) 

190  FORMATC ’ O' ,3X,’ DEMAND -BASED  PERFORMANCE  MEASURES:') 

WR ITE ( 6 , 15  0 ) SMA2 , ADDA2 , ADDB02 
WITE(6jl94J 

* '  VoViV*******  ’  ) 

RETURN 

it 

*VARIABLES  NOT  DEFINED  IN  PREVIOUS  APPENDICES 
* 

*VARIABLE  NAME  DATA 

*C0GG1  IH 

*00  REORDER  QUANTITY 

* 

END 
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APPENDIX  G.  SUBROUTINES  -  PRTOUT,  OBJECT  AND  EBO 

*THE  FOLLOWING  THREE  SUBPROGRAMS  WERE  LISTED  TOGETHER  BECAUSE 
*FUNCTION  TOGETHER 

C  ---  ROUTINE  TO  PRINT  OUT  RESULTS 

SUBROUTINE  PRTOUT(MD,NAME,QM,B,Q,QRR,QR,NT,N,NN,X,Z,Cll,D,MSRTG, 
*COG 1 , C0G2 , NPO , TOV , QMR , PBP , PVAR , SLD , PCLT , NNN , NN 1 ) 

INTEGER  ROP( 1500) ,NPO,NT,KQ( 1500) ,KQR( 1500) ,KQW( 1500) ,X( 1500) 
INTEGER  PBP(N) 

REAL  Z(N) ,0V(2,1500) ,T0V(2) ,C11(N) ,MSRTG,BSW( 1500) , SLD 
REAL  Q(N) ,D(N) ,QRR(N) ,QR(N) ,PVAR(N) 

REAL  PCLT(N) 

REAL*8  B 

CHARACTER*23  NAME 

CHARACTER*4  COGl(N) ,C0G2(N) ,QM,QMR 
COMMON  SN( 1500,9) 

DO  1  1=1, N 
KQ(I)=IFIX(Q(I)) 

ROP(I)=X(I)-KQ(I) 

BSW(I)  =  C11(I)*X(I) 

1  CONTINUE 

CALL  OBJECTCX, N,NN,Z,D,QR, PBP, PVAR, OV, TOV) 

WRITE(6j900) 

* '  *****Vf*********  '  ) 

WRITE(6,901)MD,NAME,COG1(1),COG2(1),MSRTG,QM,N,NT,NNN 

901  FORMATC  ^  0 ' , IX , ' M0DEL( ' , 1 1 , ’ ) '  IX , A2 3 , IX , ’ COG:  ' , 2A2 , 3X , 

*'MSRT  GOAL:  ' ,F8.  2,'  DAYS'  3X, ‘qP:  ' ,A4,3X, 'NI/N:  ' ,14, 

*'/' ,14, 3X, 'LOWER  BD  CODE:  V, 12) 

IF(NPO. EQ. l)GO  TO  907 
WRITE(6,902) 

902  FORMAT( 'O' ,5X, 'NIIN' ,4X, 'DEPTH' ,3X, 'MSRT(DAYS) ' ,2X, 

^'INVEST.  LVL.  ' ,5X,'C11($)' ,2X,'Q  ',3X,'ROP', 

*5X, 'PPV' ,8X, 'd' ,5X, 'PCLT' , 

*5X,'PVAR' ) 

WRITE(6,903)((SN(I,J),J=1,9),X(I),0V(2,I),BSW(I),C11(I),KQ(I), 

*ROP(I),Z(I),D(I),PCLT(I),PVAR(I),I=l,N) 

903  FORMAT(3X,9A1,2X,I5,2X,F10.4,1X,F13.  0, IX, Fll.  2,15,15, 

*F10. 3,F8.  3,1X,2X,F5.  2,F9.  1) 

906  WRITE(6,904)TOV(2),B,T0V(l),SLD 

904  FORMATC ' O' ,2X,' TOTAL  PERFORMANCE:  ' ,F9. 3, IX, '$' ,F12.  0,5X, 

*'SMA:  '  ,F8.  2, 4X,' SAFETY  LEVEL  DAYS:  ^F8.  2) 

GO  TO  909 

907  WRITE(6,908)TOV(2),B,TOV(1),SLD 

908  FORMATC 'O' ,4X,' OVERALL  PERFORMANCE:  ' ,2X,F9. 3, 4X, ,F12.  2, 5X, 
*'SMA:  ' ,F8.  2, 4X,' SAFETY  LEVEL  DAYS:  ' ,F8. 2) 

909  RETURN 
END 

* 

*  ---  ROUTINE  TO  COMPUTE  THE  OBJECTIVE  FUNCTIONS  FOR  GIVEN  ALLOCATION 
SUBROUTINE  OB JECT( X , N , NN , Z , D , QR , PBP , PVAR , OV , TOV) 
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INTEGER  N,NN,::(N),XI,PBP(N) 

REAL  Z(N) ,OV(NN,N) ,TOV(NN) ,SLT,MSRT,0V1,D(N) ,QR(N) 

REAL  CD,P,MSRTC,BO,BOT,DD,PVAR(N) 

TOV(1)=0. 

TOV(2)=0. 

B0=0.  0 
B0T=0.  0 
5  CONTINUE 
SLT=0. 

DO  10  1=1, N 

IF(Z(I).EQ.  0.  0)GO  TO  16 
XI=X(I) 

SLT  =  SLT  +  D(I) 

OV(1,I)=0, 

CALL  EBO( Z( I) , XI , D( I ) , QR( I ) , OVl , PBP( I ) , PVAR( I ) ) 
MSRTC=TWUS( Z( I ) , QR( I) ,XI , PBP( I) , PVAR( I ) .MARK) /D( I ) 
MSRT=AMAX1(MSRTC,0.  0) 

0V(1,I)=AMIN1(0V1,1. 0) 

GO  TO  18 

16  0V(1,I)=0. 0 

MSRT=1.  0 

18  0V(2,I)  =  91.*MSRT 

T0V(2)  =  T0V(2)  +  0V(2,I)*D(I) 

B0=B0+0V(1,I) 

10  CONTINUE 

IF(SLT.  EQ.  0.  0)SLT=1.  0 
BOT=BO/SLT 

T0V(1)=  (1. -B0T)*100. 

TOV(2)=TOV(2)/SLT 

RETURN 

END 


SUBROUTINE  EBO( Z , X , D , QR , OV 1 , PBP , PVAR ) 

REAL  Z,D,QR,0V1,ALPHA1,ALPHA2,X1,X2,D1,D2,P1,P2,PVAR 

REAL*8  T 1 , T2 , PHI 1 , PHI 2 , CDD 1 , CDD2 , DD 1 , DD2 , DLPHAl , DLPHA2 , PV , CD 1 , CD2 

INTEGER  PBP,X 

K1=X-IFIX(QR+.  05) 

X1=FL0AT(K1) 

X2=FL0AT(X) 

K2=X 

PV=PVAR 

IF(Z.  GT.FL0AT(PBP))G0  TO  10 
CALL  CDFP(Z,K1,P1,CD1,D1) 

CALL  CDFP(Z,K2,P2,CD2,D2) 

ALPHA1=D1*(Z-X1)+X1*P1 
ALPHA2=D2*( Z-X2)+X2*P2 
IFCALPHAl.  LT.  0.  0)ALPHA1=0.  0 
IF(ALPHA2.  LT.  0,  0)ALPHA2=0.  0 
OV 1=( ALPHA 1 - ALPHA2 ) *D/ 1 
GO  TO  20 

10  T1=(X1-Z)/DSQRT(PV) 

IF(DABS(T1).GT.  7.0)G0  TO  12 
PHIl=(DEXP(-(Tl**2)/2.  ))/SQRT(2.*3.  14159265) 

CALL  MDN0RD(T1,CDD1) 

DD1=1.  0-CDDl 
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GO  TO  13 

12  IF(T1.  GT.  7,  0)DD1=0.  0 
IF(T1.  LT.  (-7.0))DD1=1.  0 
PHI 1=0. 0 

13  T2=(X2-Z)/DSQRT(PV) 

IF(DABS(T2).GT.  7.0)GO  TO  14 
PHI2=(DEXP(-(T2**2)/2.  ))/SQRT(2.  *3,  14159265) 

CALL  MDN0RD(T2,CDD2) 

DD2=1.  0-CDD2 
GO  TO  15 

14  IF(T2.GT.  7.0)DD2=0.  0 
IF(T2.  LT.  (-7.0))DD2=1.  0 
PHI2=0. 0 

15  DLPHA 1=( PHI 1 -T1*DD 1)*( DSQRT( PV) ) 

DLPHA2=( PHI2-T2*DD2)*(DSQRT(PV) ) 

IFCDLPHAl.  LT.  0.  0)DLPHA1=0.  0 
IF(DLPHA2.  LT.  0.  0)DLPHA2=0.  0 

OV 1=( DLPHA 1 -DLPHA2 ) *D/QR 
20  RETURN 
END 

* 

*VARIABLES  NOT  DEFINED  IN  PREVIOUS  APPENDICES 
* 

*VARIABLE  NAME  DATA 
*MD  1 

*NAME  UICP  CONSUMABLES  MODEL 

*QM  (NOT  USED) 

*QRR  (NOT  USED) 

*NT  1500  (LIMIT  ON  THE  NUMBER  OF  ITEMS  THE  PROGRAM  WILL 

*  PROCESS 

*NN  2  (NUMBER  OF  MEASURES  OF  PERFORMANCE  -  SMA  &  MSRT) 

*NPO  1  OR  0  -  PRINTOUT  OPTION  (1=SUMMARY,  0=FULL  PRINTOUT) 

*TOV  COMPUTED  SMA  AND  MSRT 

*QMR  (NOT  USED) 

*SLD  DAYS  OF  SAFETY  LEVEL  FOR  THE  ENTIRE  FOUR  DIGIT  COG 

*0V  MSRT  FOR  EACH  NUN 

*0V1  EXPECTED  NUMBER  OF  BACKORDERS 

* 
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APPENDIX  H.  SUBROUTINE  CDFP 

SUBROUTINE  CDFP(ZZ,K,P,C,D) 

---  ROUTINE  TO  CALCULATE  POISSON  CDF  AND  MASS 
---  CUMULATIVE  POISSON  DISTRIBUTION 
REAL*8  ZZZ.PP,CC,CC1,DD,C 
REAL  ZZ,P,D 
INTEGER  K,I 
IF(K.  LT.  0)GO  TO  12 
ZZZ=ZZ 

PP=DEXP( -ZZZ) 

CC=PP 
CC1=0. 0 

IFCK.EQ.  0)  GO  TO  11 
KK=5*IFIX(ZZ+0. 5) 

IF(ZZ.  GT.  10.  0.  AND.K.GT.KK)GO  TO  15 
DO  10  1=1, K 
CC1=CC 

PP=PP*ZZZ/DBLE(I) 

CC=CC+PP 

10  CONTINUE 

11  P=PP 
C=CC 

DD=1.  0-CCl 
D=DD 
RETURN 
12  P=0. 0 

C=0.  0 
D=l.  0 
RETURN 
15  P=0.  0 

C=l.  0 
D=0.  0 
RETURN 
END 
C 

SUBROUTINE  CDFB(ZZ,K,PVAR,C,NB) 

C  ---  ROUTINE  TO  CALCULATE  NEGBINOMIAL  CDF 
REAL  ZZ,C,PVAR 

REAL*8  ZZZ,PP,CC,BR,R,BK,S22,B,BQ 

INTEGER  K,I,NB 

NB=0 

ZZZ=ZZ 

S22=PVAR 

BR=ZZZ/S22 

BQ=S22/ZZZ 

IF(BQ.  LE.  1.  0)G0  TO  8 

R  =  1.  0-BR 

BK=(ZZZ**2)/(S22-ZZZ) 

IF(BK*DLOG(BQ).GT.9.0)GO  TO  8 

PP=BR**BK 

CC=PP 

IF(K  .EQ.  0)  GO  TO  11 
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GO  TO  9 

8  NB=1 
RETURN 

9  DO  10  1=1, K 
B=DBLE(I-1) 

PP=PP*R*( B+BK) /DBLEC I ) 

CC=CC+PP 

10  CONTINUE 

11  C=CC 
RETURN 
END 

* 

*ALL  VARIABLES  EXCEPT  OUTPUT  VARIABLES  (P,  NB)  PREVIOUSLY  DEFINED 
* 
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APPENDIX  I.  SAMPLE  CHURN 
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